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ABSTRACT 

We trace the tidal Stream of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using Red Clump 
stars from the catalog of the Sloan Digital Sky Survey- Data Release 6, in the range 150° < RA < 220°, 
corresponding to the range of orbital azimuth 220° < A < 290° . Substructures along the line of sight 
are identified as significant peaks in the differential star count profiles (SCP) of candidate Red Clump 
stars. A proper modeling of the SCPs allows us to obtain: (a) < 10% accurate, purely differential 
distances with respect to the main body of Sgr, (b) estimates of the FWHM along the line of sight, 
and (c) estimates of the local density, for each detected substructure. In the range 255° ;^ A < 290° 
we cleanly and continuously trace various coherent structures that can be ascribed to the Stream, in 
particular: the well known northern portion of the leading arm, running from d ~ 43 kpc at A ~ 290° 
to d ~ 30 kpc at A ~ 255°, and a more nearby coherent series of detections lying at constant distance 
d ~ 25 kpc, that can be identified with a wrap of the trailing arm. The latter structure, predicted 
by several models of the disruption of Sgr dSph, was never traced before; comparison with existing 
models indicates that the difference in distance between these portions of the leading and trailing arms 
may provide a powerful tool to discriminate between theoretical models assuming different shapes of 
the Galactic potential. A further, more distant wrap in the same portion of the sky is detected only 
along a couple of lines of sight. For A < 255° the detected structures are more complex and less 
easily interpreted. We are confident to be able to trace the continuation of the leading arm down 
to A ~ 220° and d ~ 20 kpc; the trailing arm is seen up to A ~ 240° where it is replaced by more 
distant structures. Possible detections of more nearby wraps and of the Virgo Stellar Stream are also 
discussed. These measured properties provide a coherent set of observational constraints for the next 
generation of theoretical models of the disruption of Sgr. 

Subject headings: galaxies: dwarf- Galaxy: structure - (galaxies:) Local Group - stars: distances — 
Galaxy: formation 



1. INTRODUCTION 

Stellar tidal streams as well as other substructures in 
the Milky Way (MW) halo are generally interpreted as 
the relics of the process of hierarchical formation of the 
MW, as envisaged by the currently accepted cosmolog- 
ical model (A-Cold Dark Matter, A-CDM hereafter, see 
Bullock et al. 2001; Madau et al. 2008, and references 
therein). With the advent of large modern surveys, like 
the 2 Micron All Sky Survey (2MASS, Skrutskie et al. 
2006) and the Sloan Digital Sky Survey (SDSS, Adelman- 
McCarthy et al. 2008, and references therein), our abil- 
ity to detect stellar systems and/or structures in the 
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halo and in the disk of the MW has increased dramati- 
cally and several large-scale likely relics of the build-up 
of the Galactic halo have been identified (Ibata et al. 
2001a; Newberg et al. 2002; Yanny et al. 2003; Majew- 
ski et al. 2003; Martin et al. 2004; Belokurov et al. 2006; 
Juric et al. 2008). Also smaller tidal streams have been 
found around disrupting globular clusters (see, for ex- 
ample Rockosi et al. 2002; Grillmair & Johnson 2006) 
or lacking an evident progenitor (Grillmair & Dionatos 
2006; Belokurov et al. 2006, hereafter Bel06). The most 
spectacular example of the process of tidal disruption and 
accretion of a dwarf satellite into our Galaxy is the Sagit- 
tarius dwarf spheroidal galaxy (Sgr dSph), originally dis- 
covered by Ibata et al. (1994). The main body of the Sgr 
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galaxy is located at ^ 26 kpc (Monaco et al. 2004) from 
the Sun, beyond the Galactic bulge (Galactic coordinates 
;, 6 = +5.6°, -14.0°). The stellar content of the Sgr dSph 
is dominated by an intermediate-age relatively metal-rich 
population, with distributions peaking at age ^ 6-8 Gyr 
and [Fe/H] ~ —0.5, (see Bellazzini et al. 2006a, here- 
after B06a, and references therein) but there is also clear 
evidence for the presence of an older (> 10 Gyr) and 
more metal-poor population as well, including Blue Hor- 
izontal Branch (BHB, Ibata et al. 1997; Bellazzini et al. 
1999a; Monaco et al. 2003) and RR Lyrae stars (Mateo 
et al. 1995a; Alcock et al. 1997; Gseresnjes 2001). All the 
available spectroscopic analyses indicate that the Metal- 
licity Distribution (MD) of Sgr stars is characterized by 
a broad peak in the range —1.0 < [Fe/H] < 0.0, with 
a weak tail likely extending beyond [Fe/H] < —2.0 (see 
B06a, Gseresnjes 2001; Monaco et al. 2005a; Mc William 
& Smecker-Hane 2006; Bonifacio et al. 2006; Sbordone 
et al. 2007; Bellazzini et al. 2008; Lagadec et al. 2009). 

The body of Sgr dSph appears tidally disturbed (Ibata 
et al. 1995), and, soon after its discovery, it was real- 
ized that there was some tidal debris surrounding the 
galaxy (Mateo et al. 1996; Fahlman et al. 1996; Alard 
1996; Ibata et al. 1997; Majewski et al. 1999). Indeed, 
it has been subsequently established that there are two 
huge tidal tails emanating from the edges of the galaxy 
and approximately tracing its orbital path, as expected 
from N-body simulations (Johnston et al. 1995; Ibata & 
Lewis 1998). These tails form a coherent and dynami- 
cally cold filamentary structure (hereafter Sgr Stream) 
that extends for tens of kpc from the parent galaxy and 
has been probed with many different tracers. Yanny 
et al. (2000) used SDSS first-year commissioning data 
to identify an overdensity of blue A-type stars in two 
stripes located at {I, b, D)= (341°, -^57°,46 kpc) and 
(157°, —58°, 33 kpc), which were subsequently found to 
match the prediction of the available theoretical models 
of the Sgr Stream (Ibata et al. 2001a). Similarly, Ivezic 
et al. (2000b) noticed that clumps of RR Lyrae stars in 
SDSS commissioning data lay along the Sgr orbit. The 
thorough study of the structure of the halo as traced by 
F stars from the SDSS, within a strip of ±1.26° around 
the celestial equator, by Newberg et al. (2002) provided 
the first examples of Golor Magnitude Diagrams (GMD) 
of the Stream population toward {I, b)= (350°, 50°) and 
(/, h)= (167°, —54°). Other detections toward specific di- 
rections were provided by Martinez-Delgado et al. (2001, 
2004), Bellazzini et al. (2003b) and Vivas & Zinn (2006). 
The first panoramic view of the Sgr Stream was presented 
by Ibata et al. (2002) using late M giants selected from 
the Second Incremental Data Release of 2MASS. Subse- 
quently, Majewski et al. (2003, hereafter M03), having at 
disposal the final All Sky Data Release of 2MASS, used a 
larger sample of M giants to provide a clearer view of the 
whole complex, tracing very neatly the trailing tidal tail 
all over the Southern Galactic hemisphere, as well as the 
part of the leading arm closer to the main body of the 
galaxy, up to RA ~ 190°. In a very recent analysis Yanny 
et al. (2009a) showed that M and K giants can be success- 
fully selected also from the SDSS and SEGUE datasets 
(Yanny et al. 2009b) and used to trace the Stream; one 
main advantage of using giants as tracers is that they can 
be (relatively) easily foUowed-up spectroscopically, thus 
providing crucial kincmatical and chemical information 



(Majewski et al. 2004; Monaco et al. 2007; Ghou et al. 
2007; Yanny et al. 2009a). 

Bel06 exploited the SDSS data release 5 (DR5) to pro- 
vide a picture of the leading arm of the Sgr Stream in the 
vicinity of the North Galactic Gap with remarkable clar- 
ity, using tracers (A-F dwarfs) that are intrinsically more 
numerous than M giants, for a given space density and/or 
surface brightness ^ . In their Fig. 1 they show the density 
of (candidate) A-F dwarf stars (selected with a simple 
color cut, g ~ r < 0.4, corresponding to ^ B ~V < 0.6) 
in the portion of the sky covered by the SDSS. The Sgr 
Stream emerges very clearly as a broad (and bifurcated) 
band going from (a,5)~(220°, 0°) to (a,5)~(125°, 25°), 
where it plunges into the Galactic Disk. The color cut 
adopted by Bel06 is very successful in tracing the Stream 
structure as it takes advantage of the fact that Sgr stars 
are younger than typical halo stars, hence they have a 
bluer Turn Off (TO) color with respect to the halo pop- 
ulation (see Unavane et al. 1996, for another application 
of the same principle). The density map by Bel06 shows 
evidence for a clear distance gradient along the Stream, 
from the nearest part crossing the Disk at a « 120°, to 
the most distant part at a « 210°, toward the North 
Galactic Pole (NGP). More recently, in a pilot project 
limited to a sub-sample of the SDSS (the so called Stripe 
82) Gole et al. (2008) described a more refined approach 
to the study of the spatial structure of the Stream, us- 
ing the same tracers as Newberg et al. (2002). Very re- 
cent detections from different data and/or using different 
models can be found also in de Jong et al. (2009), Prior 
et al. (2009b) and Keller (2009a). 

In Bellazzini et al. (2006c, hereafter B06c) we demon- 
strated that yet another kind of tracer can be efficiently 
used to study the Sgr Stream, i.e. core-He-burning stars 
lying in the well populated Red Glump (RG) of the GMD 
of Sagittarius dSph. In particular, wc showed that it is 
possible to detect the RG associated with a given sub- 
structure as a peak in the differential Star Gount Pro- 
files (SGP) of sub-samples of stars selected in a rela- 
tively narrow color range including the RG. The spatially 
localized RG population can be disentangled from the 
fore/background contaminating population of the MW 
by subtracting the underlying SGP, that is, in general, 
quite smooth and smoothly varying with position in the 
sky. In B06c wc used this technique to compare the Hor- 
izontal Branch (HB) morphology in the Stream and in 
the main body of Sgr, finding an age/metallicity gradient 
along the Sgr remnant (see also Monaco et al. 2007; Ghou 
et al. 2007), while in Gorrenti et al. (2009) we obtained 
an independent detection of the recently discovered stel- 
lar system Bootes HI (Grillmair 2009; Garhn et al. 2009), 
providing new insight on its nature, structure and stel- 
lar populations. Garrell & Wilhelm (2010) recently pre- 
sented the results of a spectroscopic survey targeting RG 
stars in the Sgr Stream, selected as in B06c. 

^ In a stellar population of given age and chemical composition 
the number of stars per units of sampled (integrated) luminosity 
in a given evolutionary phase is proportional to the duration of the 
evolutionary phase (see Renzini & Buzzoni 1986; Rcnzini & Fusi 
Pecci 1988; Renzini 1998, and references therein). A-F stars are 
evolving along the Main Sequence, a phase lasting several Gyr for 
these stars, while M giants are in the latest phases of their evolu- 
tion along the Red Giant Branch, lasting < 10* yr. Hence in any 
given field, independently of the absolute density normalization, 
A-F dwarfs outnumber M giants by a factor of > 10. 
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The most natural and direct application of this tech- 
nique is the determination of accurate distance estimates 
from the magnitude of detected RC peaks, as the RC 
is well known and widely used as a standard candle 
since long time (see Paczynski & Stanek 1998; Stanek 
& Garnavich 1998; Girardi & Salaris 2001; Babusiaux 
& Gilmore 2005; Bellazzini et al. 2006b, and references 
therein). For intermediate/old-age populations, the lu- 
minosity of the RC peak shows relatively modest vari- 
ations as a function of age and metallicity, in particu- 
lar when measured in the reddest optical passband (as 
Cousins'/, see Girardi & Salaris 2001). When used dif- 
ferentially, i.e. looking at the same (or very similar) stel- 
lar population in different places, the variations in the 
intrinsic luminosity of the RC due to age/metallicity ef- 
fects should vanish. Given also the intrinsic narrowness 
of the feature in Sgr (see below) , the RC seems the ideal 
tool to accurately trace the run of the distance along the 
orbital azimuth of the Sgr Stream, from the main body 
of the galaxy all over the portion of the Stream sampled 
by the SDSS. The other large survey covering all the ex- 
tent of the Stream, 2MASS, cannot be used in this way 
as the associated photometry is not sufficiently deep to 
reach the RC level. 

In this paper we will use the RC SCP method outlined 
above to take accurate purely differential measures of the 
distance of the Northern arms of the Sgr Stream with re- 
spect to the main body of the galaxy. This will provide 
strong constraints for the models of the disruption of Sgr 
within the Galactic (dark) halo, and, in turn on the phys- 
ical properties of the dark halo itself (Ibata et al. 2001b; 
Hclmi 2004; Johnston et al. 2005). The basic idea is the 
following: (a) we measure the position of the RC peak in 
V and /, with independent color selections using B — V 
and V — I colors, in the main body of Sgr (from B06a 
photometry), (b) we select SDSS fields projected onto 
the Sgr Stream as traced by Bel06, (c) we transform the 
SDSS photometry into B,V,I magnitudes, (d) we de- 
tect the RC peak(s) in V and / SCPs from the SDSS 
on-Stream fields (adopting the same color selections as 
in the main body), and (e) we obtain two independent 
measures of the magnitude differences of the RC peaks 
between the main body and the considered portion of 
the Stream. These are fully equivalent to differences in 
distance modulus, that is, differences in distance. The 
whole set of differential distances can be translated into 
a set of absolute distances by adopting the preferred value 
of the distance modulus for the main body (see, for ex- 
ample Alard 1996; Layden & Sarajedini 2000; Monaco 
et al. 2004; Kunder & Chaboyer 2009). The detection of 
the same peaks in both V and / SCPs provides a useful 
sanity check on the interpretation of the SCPs and on 
the derived differential distances. As an additional ob- 
servational constraint to models of the disruption of Sgr, 
we provide also an estimate of the characteristic width of 
the Stream section crossed by our fields (see Sect. 3.5). 

The plan of the paper is the following. In Sect. 2 we 
present the field of the main body that we used as tem- 
plate and the fields of the Stream used for/in the anal- 
ysis. In Sect. 3 we describe the method used to analyze 
the SCPs and derive the informations from the peaks. In 
Sect. 4 we present all the SCPs obtained from each field 
and we discuss some special cases. In Sect. 5 we compare 
our results with previous works in literature, with par- 



ticular emphasis on the different degree of uncertainty 
related to the distance estimates. In Sect. 6 we com- 
pare our distance estimates with models that reproduce 
the three-dimensional shape of the Stream. Finally, we 
summarize and discuss our results in Sect. 7. 

Some preliminary reports on earlier phases of this 
project were presented in Correnti et al. (2007) and Cor- 
renti et al. (2008). 

2. DATA AND OBSERVABLES 

As a reference sample for the stellar population in the 
core of Sgr we take the photometry of a 1° x 1° wide 
field located ^ 2° eastward of the galaxy center at {l,b) 
~ (6.5°, —16.5°), presented in B06a and named Sgr34- 
This should be considered as fairly representative of 
the average population of the Sgr galaxy (see Bellazzini 
et al. 1999a, b; Giuffrida et al. 2010), avoiding the 
youngest and most metal-rich populations that appear 
to reside in the central nucleus (Siegel et al. 2007; 
Bellazzini et al. 2008). The strong similarity between 
the population of the Sgr main body and the Stream 
has been shown by Newberg et al. (2002) and Bcl06, by 
direct comparison of CMDs^. To sample the Galactic 
population at similar angular distance from the Galactic 
Center as for Sgr34 we used the same control field also 
presented in B06a: a 0.5° x 0.5° field, named GaLField, 
at {l,b) ~ (-6.0°, -14.5°), that was used in B06a to 
perform the statistical decontamination of the Sgr34 
CMD from the foreground/background Galactic stars. 
Following B06a we adopted the average reddening values 
{E{B - V)) ^ 0.116 for Sgr34 and {E{B - V)) ^ 0.096 
for GaLField, as derived from the reddening maps of 
Schlegel et al. (1998, hereafter SFD98). 

To study the Stream, we used the SDSS-DR6 photome- 
try of objects classified as stars (extracted from the SDSS 
CasJobs query system, Adelman-McCarthy et al. 2008)"^ 
for a series of selected fields along the branch A and B, 
listed in Tab. 1 and plotted in Fig. 1. We chose to follow 
the two branches separately, with non-overlapping fields. 
These on-Stream fields are similar, in position, to those 
studied by Bel06, but are slightly smaller (5° x 5° instead 
of 6° X 6°), to avoid overlap between different fields of the 
same branch. For each on-Stream field [located, for ex- 
ample, at (l,b)=(^Oj bo)] we selected also a corresponding 
control field (CF) located at the same latitude and at the 
same angular distance from the Galactic Center on the 
other side of the Galaxy [i.e having (l,b)=(360° — ^o, ^o)]^- 
Assuming that the MW is symmetric about its center 
and its disk mid-plane (that should be a reasonable first- 
order approximation, at least at the Galactic latitudes 
considered here, b > 45.5°; but see Bell et al. 2008), each 
control field should be fairly representative of the typ- 
ical Galactic population contaminating our on-Stream 
fields. Following Bel06, to average out the effects of shot 
noise, the control fields are larger than the on-Stream 

^ In particular, Bcl06 uses the same photometry of Sgr34 that 
is adopted here, as a reference. 

^ We used a template SQL query provided in the Cas- 
SDSS web page called "Clean Photometry", in the ver- 
sion aimed at selecting stars. The corresponding SQL 
string can be found under the link Clean Photometry at 
http : //cas . sdss . org/astrodr7/en/help/docs/real query . asp. 

"^ Except in the case of field FIB, for which the CF is the same 
adopted for the field F2B 
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Fig. 1. — Left panel: Distribution of the on-Stream fields of the branch A (red squares) and the braneh B (blue squares) of the Sgr 
Stream. The associated control fields are also plotted (black dashed squares). The plot is intended to show the position of the fields and 
give a rough idea of their dimension, projection effects are not taken into account. The shape of the two branches is reproduced with a 
continuous line, following Fig. 1 of Bel06. The positions of known stellar systems falling into the considered fields are indicated (green 
stars) and labeled. Right panel: Positions of the various fields in Galactic coordinates, to highlight the symmetry (with respect to the 
Galactic Center and to the Galactic Plane) of the adopted on-Stream fields and their corresponding Control Fields. 

TABLE 1 
Position and reddening of the considered fields'*. 
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51 


70 
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40 


S3 
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62 
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62 
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17 


224.5 
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0.030 
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60 
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0.014 
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14A 
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50 
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0.004 
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62 
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50 
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0.006 


15A 
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19 
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45.5 


0.030 


0.006 
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63.5 
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45.5 


0.017 
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IB 
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3B 


205 


9 
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6B 
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81 


0.025 
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33 


75 


81 


0.012 


0.002 


7B 


185 


20 
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80 


0.029 


0.005 
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80 


0.013 


0.003 


8B 
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22 
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77 


0.028 


0.006 
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40 


126 


77 


0.016 


0.004 


9B 


175 


22.75 


224 


73 


0.023 


0.004 


187.5 


43 


136 


73 


0.017 


0.005 


lOB 


170 


24 


216.5 


69 


0.017 


0.002 


182.5 


46.5 


143.5 


69 


0.017 


0.005 
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165 


25.5 


210.5 


65 


0.022 


0.008 


176 


48.5 


149.5 


65 


0.019 


0.005 


12B 


160 


25.5 


208.5 


60.5 


0.025 


0.007 


175 


48 


151.5 


65 
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0.005 


13B 


155 


27.5 


203 


56.5 


0.028 


0.007 


163 


52 


157 


56.5 


0.013 
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^{E{B — V)) is the mean reddening of the field as extracted from 
Schlegel et al. (1998) maps and averaged over all the stars in the 
field; a is the corresponding standard deviation. The subscript c 
refers to control fields. 

fields (10° X 10°). As shown in Fig. 1, the globular clus- 
ters NGC5466 and M3, and the dwarf galaxy remnant 
Bootes III are enclosed within some of our control fields 
(Correnti et al. 2009). To avoid any undesired contami- 
nation we excluded from the corresponding samples the 
stars associated with these stellar systems by excising ar- 
eas of radius 1° (for the globulars) and 2° (for the dwarf 
galaxy) around their centers. The only known stellar sys- 
tem that is (partially) enclosed in one of our on-Stream 
fields is the globular cluster NGC4147 (see Bel06 and 
Bellazzini et al. 2003a, b). Also in this case we excluded 
from the adopted sample all the stars within 1° of the 
cluster center. In the following we will use CFs only as 
a further observational check that the simple models we 



adopt to account for the fore/background populations 
contaminating our SCPs are adequate for our purposes 
(see Sect. 3.1). 

For our analysis we adopted the reddening-corrected 
g, r, i, z magnitudes as provided by CasJobs. These mag- 
nitudes were also corrected using the SFD98 maps, hence 
the source of the reddening corrections is homogeneous 
for all the datasets considered in the present analysis. 
The mean E{B — V) and its standard deviation for each 
field, averaged over all the stars included in the field, 
are reported in Table 1. It is important to note that 
the average reddening of our fields is remarkably low 
(0.010 < E{B -V) < 0.038) and constant within each 
field (0.002 < (7e{b-v) < 0.010), hence any error in 
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the adopted reddening correction would have only a mi- 
nor impact on our final differential distance estimates. 
For brevity, in the following all the reported magni- 
tudes and colors are reddening-corrected (i.e., for exam- 
ple, V — Vq ~ extinction corrected V magnitude). 

The g, r, i, z magnitudes in the SDSS system have been 
transformed to the Johnson-Kron-Cousins B, V, I system 
(as defined by the standard stars by Landolt 1992) using 
robust empirical transformations that have been checked 
to be particularly accurate in the color range typical of 
RC stars (provided by Lupton 2005^, derived from large 
samples of stars in common between SDSS and the ex- 
tended database of Landolt's standards by Stetson 2000). 
In particular, we obtain B and V from g and r, while / is 
obtained from i and z, adopting the following equations: 



B = 


= g + 0.3130(.g - r) + 0.2271 


(J - 


= 0.0107 


(1) 


V = 


= g - 0.5784(5 - r) - 0.0038 


a - 


= 0.0054 


(2) 


I = 


= i-0.3780(i-z) -0.3974 


a - 


= 0.0063 


(3) 



Note that the transformed V and / are fully inde- 
pendent as they are obtained by independent couples of 
SDSS magnitudes. Consequently, measures of the po- 
sition of any significant peak detected in V and / SCPs 
will also be independent, thus providing a powerful cross- 
check of any detection and distance estimate. 

2.1. Selections on the Color Magnitude Diagram 

In Fig. 2 we present reddening-corrected V,B-V and 
I, V-I CMDs (focused on the RC features up to the upper 
region of the Main Sequence, MS) of the main-body field 
Sgr34 and of the corresponding control field GaLField. 
The comparison between the CMDs of the two fields 
permits the identification of the main features associ- 
ated with Sgr and with the fore/background Galactic 
populations. The RC of the Sagittarius dSph is a promi- 
nent feature in the CMDs of the Sgr34 field (upper pan- 
els), around {I,V - I) ~ (16.9,0.9) and iV,B - V) ~ 
(17.8,0.8). The wide and inclined Red Giant Branch 
(RGB) can be discerned over the background, going 
from (/, V - I) c^ (16.9, 0.9) to (/, V - I) ^ (14.0, 1.5) 
[iV,B - F) ~ (19.5,0.8) to {V,B - V) ~ (16.0,1.4)], 
and continuing beyond the limits of the box. The RGB 
bump is apparent at T^ ^ 18.2 and / ^ 17.2, along 
the RGB (see Monaco et al. 2002). For y - / < 0.3 
{B -V < 0.2) at / - 19.0 {V - 17.9) a portion of the 
Blue Horizontal Branch (BHB) is also visible (Monaco 
et al. 2003); a.t V - I <0.5 {B -V < 0.4) and / > 19.5 
(y > 18.5) the Blue Plume (BP, Mateo et al. 1995b; Bel- 
lazzini et al. 1999a, B06a) population is visible. The Sub 
Giant Branch (SGB, for ^ - / > 0.8 or B - V > 0.7) and 
the upper Main Sequence (MS, to the blue of the above 
limits) appear for I > 19.0 {V > 20.0). For a more 
detailed description of the CMD of Sgr see B06a. The 
strong vertical band around V — I ^ 0.7 {B ~ V ^ 0.6) 
running over the largest part of the CMD, and bending to 
the red at / ~ 19, F ~ 20, is constituted by MS stars of 
the MW (mostly from the Thick Disk, in this direction, 
according to the Galactic model by Robin et al. 2003); 
the wide band running parallel, to the red of the vertical 
portion of this feature is mainly populated by Galactic 

^ http://www.sdss.org/dr7/algorithrns/sdssUBVRITransfomi.html 



giants, either in their RGB or RC/HB phase. The ma- 
jority of the stars redder than V~I = B — V^ 1.0 
belongs to the vertical plume of local Galactic M dwarfs. 

The vertical lines in each panel of Fig. 2 enclose the 
color stripes that we adopted to select the RC population 
in the two colors, corresponding to 0.70 < B ~V < 0.95 
and 0.85 < V — I < 1.05. The choice of the color lim- 
its was made in order to include the bulk of the RC 
population even if small color shifts were present due to 
errors in the adopted reddening corrections and/or pop- 
ulation gradients, while keeping the contamination from 
other sources as low as possible. The distribution in color 
within the selection windows (around the magnitude of 
the observed RC of Sgr) shown in Fig. 3 suggests that 
color shifts of order ±0.05 mag would lead just to mi- 
nor losses of the signal (of the order of 10% with respect 
to the number of stars obtained with our choice in the 
selection window) . 

Fig. 2 clearly shows that, in addition to Sgr RC stars, 
several different contaminants are expected to enter the 
selection window. For / < 18.5 {V < 19.5) Galactic gi- 
ants (mainly RC stars) should be the primary source of 
contamination, while the sequence of Galactic MS stars 
crosses the selection stripes at / > 19 and V > 20, boost- 
ing the star counts at faint magnitudes. The RGB of the 
Sgr population, and in particular the RGB bump, are 
also selected by the adopted windows. We will show be- 
low that this source of contamination has a negligible 
effect on our SCPs. At / > 19.0 {V > 20.0) the SGB 
stars of Sgr enter the selection window; as they are much 
more numerous than RGB and RC stars they may pro- 
vide a serious contribution to the "background" in our 
SCPs, at faint magnitudes. Finally, the MS of Sgr crosses 
the windows at / > 21 (T^ > 22). The actual structure 
of the contamination entering the windows will obviously 
depend on (a) the Galactic population encountered along 
the considered line of sight {los , hereafter), and (b) the 
distance of the wrap(s) of the Sgr Stream that is(are) 
crossed by the considered los . However the los consid- 
ered in this study are all at much higher Galactic lat- 
itudes than Sgr34, hence the degree of contamination 
per unit area of the sky should be lower, and the average 
distance of the encountered stars should be higher, hence 
most of the contamination by Galactic dwarfs should oc- 
cur at fainter magnitudes than discussed above for Sgr34- 
Furthermore, all the detections of the Stream presented 
here are at distances similar or larger than the main body 
of the galaxy sampled by Sgr 34', hence, in most cases the 
contamination by the SGB of the Stream population will 
occur at fainter magnitudes than in Sgr34^ ■ In any case, 
to limit the contribution by dwarf stars, independently 
of their origin, we limit our analysis to the magnitude 
ranges 15.0 < I < 19.5 and 16.0 <V< 20.5^. These 
limits approximately correspond to an accessible range of 
heliocentric distances 12 kpc ^ -D < 70 kpc (see Fig. 5, 
below). 

^ However it should be noted that the MS of wraps of the Streams 
that are too nearby to have their RC detected with the present tech- 
nique may contribute to the contamination of our color-selected 
samples of candidate RC stars. Moreover, other unknown sub- 
structures may contribute to the contamination (see Correnti et al. 
2009, for example). 

^ Except in the case of Sgr34 where the limits arc 15.0 < I < 18.5 
and 16.0 <V < 19.5 
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Fig. 2. — Reddening corrected {V,B-V; left panels) and {I, V-I; right panels) CMDs, focused on the RC up to the upper part of the 
MS, of the Sgr field (Sgr34, upper panels) and of the control field (GaLField, lower panels). The vertical lines enclose the RC populations 
and are the color strips used to select the region where build the SCPs. The color ranges are respectively, 0.70 < {B — V) < 0.95 and 
0.85 < (V - /) < 1.05. 



While the surface brightness of Sgr at Sgr34 is ^ 25 
mag/arcsec^, typical values for the Stream are > 30 
mag/arcsec^ (Mateo et al. 1998; Bellazzini et al. 2003b; 
Majewski et al. 2003, and references therein). It may be 
quite hard to identify the feeble signal from such sparse 
populations even in the presence of low background. In 
fact, even in the most favorable cases, the RC is barely 
visible in the CMDs of on-Stream Fields (see, e.g. New- 
berg et al. 2002) . The construction and modeling of SCPs 
described below is very effective in extracting the dis- 
tance information in these cases (B06c, Correnti et al. 
2009). 

Finally, there are several indications that there is a 
sizable metallicity (and presumably age) gradient along 
the Stream, in the sense that the average metallicity is 
lower in distant portion of the Stream with respect to 
the main body of Sgr (B06c, Monaco et al. 2007; Chou 
et al. 2007). This is generally interpreted as due to a pre- 
existing population gradient within the original body of 
the Sgr galaxy, as the tidal tails were preferentially pop- 
ulated by stars that resided in the old and metal poor 
outskirts of Sgr (Chou et al. 2007). It must be stressed 
that the detected gradient means that the relative pro- 



portion of intermediate-age & metal-rich stars and of old- 
age metal-poor stars changes along the Stream (and with 
respect to the main body). This, in turn, changes the HE 
morphology, i.e. the relative abundance of RC and Blue 
HB stars (as observed in B06c) , but is it not expected to 
change the intrinsic luminosity of the RC. Indeed, Car- 
rell & Wilhelm (2010) find that the mean metallicity of 
RC stars along the Stream is very similar to that found 
in the main body of Sgr. Hence, while the population 
gradient may bias estimates of the stellar density along 
the Stream obtained from RC stars, our distance esti- 
mates should be unaffected and our characteristic size 
estimates can be only marginally affected (see Sect. 3.5 
and Fig. 18, below, for further details and discussion). 

2.2. Detecting RC peaks in Star Count Profiles 

All the SCPs of color-selected RC samples presented 
in this paper are computed as running histograms (see 
Bellazzini et al. 2005, and references therein)*, as these 
couple the property of collecting the signal from a wide 

* Running histograms are histograms in which the step is smaller 
than the bin width. The adoption of steps much smaller than the 
bin width removes the dependency from the starting position of 
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Fig. 3. — Distribution of stars in B-V (upper panel) and V-I 
(lower panel) eolor stripes for the Sgr34 field. The selection must be 
optimized in order to contain all the RC stars, but minimizing the 
presence of contaminating stars. The RC is well confined inside the 
selected regions and peaks near the middle in both; the magnitude 
limits of the selection are respectively, 17.3 < V < 18.3 and 16.4 < 
/ < 17.4. 

bin with the abihty of constraining with great accuracy 
the location of density maxima, almost independently of 
the bin width. A bin width of ±0.2 mag^ and step of 
0.02 mag have been adopted here. After different tri- 
als, they have been found to provide a good compro- 
mise between the exigence of co-adding all the signal 
from a given RC population (that requires larger bins) 
and the ability of distinguishing (resolving) nearby peaks 
(that is favored by smaller bins). The use of generalized 
histograms (Laird et al. 1988) would have provided an 
higher degree of smoothing, possibly making some of our 
SCPs easier to interpret. However we preferred running 
histograms as they provide the reader a clearer idea of 
the local noise on the SCP as well as a scale in real units 

( 1*-°''^F 1/ )- The density scales of the various fields 

have all been reported to unit standard area (1° x 1°) by 
applying the corrections due to spherical geometry that 
is inherent to equatorial coordinates. 

To illustrate at best the case of the detection of the RC 
of a spatially confined stellar system in a color-selected 
SCP, we show in Fig. 4 the V and / SCPs for the Sgr34 
field (continuous lines), compared with those obtained 
for the control field GaLField, normalized by the ratio of 
background densities between the two fields^°. 

The shapes of the Sgr34 and GaLField SCPs are re- 
markably similar: the only exception is the very strong 
and well defined peak corresponding to the RC of the 
Sgr galaxy seen in Fig. 2. It is interesting to note that 

the binning that affects classic histograms. Clearly, the values of 
adjacent bins are not statistically independent. 

^ With rare exceptions in which bins of ±0.25 mag have been 
adopted to enhance the signal of a weak feature. All these cases 
are clearly indicated in the following. 

^"^ This ratio is dominated by the ratio of the areas of the fields, 
3gr34 being ~ 4 times larger than GaLField. However GaLField 
sample a direction ~ 2° closer to the Galactic Plane and ~ 0.5° 
closer to the Galactic Center than Sgr 34, hence the (column) stellar 
density is intrinsically larger in the former field. The actual density 
ratio, computed in selected CMD boxes where the contribution 
from Sgr dSph stars is negligible, is ~ 3, see Bellazzini et al. (1999b) 
and B06a. 




Fig. 4. — De-reddened running histogram SCPs of color-selected 
RC candidates for Sgr34 (continuous line) and for the control field 
GaLField (dashed line) in V and / magnitude (respectively upper 
and bottom panel). The shape of the SCPs is similar but, as ex- 
pected, in the GaLField the peak is completely laking, while in the 
Sgr field it stands out very clearly. 

while also other features related to Sgr are visible in the 
CMDs and (at least partially) included in the selection 
windows, as for example the RGB bump, in the SCPs the 
RC is the only signal emerging from the Sgr population. 
Independently of the origin of the stellar mix actually se- 
lected, the SCP of the control field, and, by analogy, the 
SCP of the contaminating population that is superposed 
on the RC in the Sgr34 SCP, are quite smooth and have 
a very simple behavior; in B06c, Correnti et al. (2009) 
and in Sect. 3.1 we show that this is the general behavior 
of the SCP of the back/foreground population in the vast 
majority of the considered los , thus justifying the choice 
of a very simple model for them, as described in Sect. 3, 
below. 

2.3. Sensitivity of the technique 

Before proceeding with the description of the method 
adopted to obtain the actual differential distance esti- 
mates, it may be useful to study the sensitivity of our 
SCPs to the various properties of any encountered sub- 
structure (distance, density, etc.). To do that we used the 
dedicated web tool^^ of the BASTI repository of stellar 
models (Pietrinferni et al. 2004; Cordier et al. 2007) to 
produce a synthetic population of ~ 45000 stars having 
age and metallicity similar to the bulk of the Sgr popu- 
lation (age=6 Gyr, [Fe/H] = -0.5). The CMD and the 
color-selected RC SCP of the population are shown in the 
left panels of Fig. 5. The synthetic stars have been dis- 
tributed along the line of sight according to gaussian dis- 
tributions having various mean distances (D=15, 25, 40 
and 65 kpc) and Full Width at Half Maximum (FWHM) 
of 3.3 kpc or 6.6 kpc, to simulate the crossing of a wrap of 
the Stream at various distances and with different char- 
acteristic sizes along the los . A FWHM ~ 3.3 kpc is 
quite typical of sections of Stream wraps crossed perpen- 
dicularly by a given los , as measured on the models of 
the disruption of the Sgr galaxy by Law et al. (2005). 
The FWHM ~ 6.6 kpc has been considered to account 

^^ http : //alblone . oa-teramo . Inaf . it 
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Fig. 5. — (a). Left panel: CMD of a synthetic population composed of 45953 stars having age=6 Gyr, [Fe/ H\ = —0.5, solar scaled 
abundance pattern, 10% of binary systems, and Salpeter's Initial Mass Function, from the BASTI database. The color window adopted 
to select RC candidates is enclosed within the two thin parallel lines. Right panel: SCP of the color-selected RC candidates obtained in 
the same way as Fig. 4(x symbols) and fitted with the same kind of model (red continuous line). The mean and standard deviation of the 
best-fitting Gaussian are reported within parentheses, (b). The synthetic population shown in panel ('aj has been distributed along the line 
of sight according to gaussian distributions having various mean distances (D=15, 25, 40 and 65 kpc, from the upper to the lower panel, 
respectively) and Full Width at Half Maximum of 3.3 kpc (thick lines) or 6.6 kpc (thin lines). The SCP of the resulting color-selected 
RC population has been derived and superposed to the (model) SCP of the background of Fig. 4, to simulate the detection of the same 
structure at different distances with the method applied here. The vertical segments indicate the positions of the peaks produced by the 
considered structures. 



for cases of sparser portions of the Stream and/or non- 
perpendicular intersections with the los . The SCPs of 
the resulting color-selected RC population have been de- 
rived (properly including realistic photometric errors^^) 
and added to the SCP of the background of Fig. 4, as 
modeled in Fig. 6, below, to simulate the detection of 
the same structure with the method applied here. The 
results of this simulation are shown in the right panels of 
Fig. 5. 

The most obvious effect shown in Fig. 5 is the increase 
of sensitivity with the distance of the structure. This is 
due to two factors: (1) the inherent "compressive" prop- 
erty of the magnitude scale, by which, for instance, a 
difference in distance of 3.3 kpc corresponds to a differ- 
ence in magnitudes of 0.43 mag at D=15 kpc and to just 
0.14 mag at D=50 kpc, and (2) the relative dimension 
of the considered structure and of the los cone at the 
distance of the structure; for nearby structures the fixed 
FoV adopted here may be smaller than the structure it- 
self, thus missing part of the signal that instead would 
be included when more distant structures are encoun- 
tered. This effect illustrates a fundamental property of 
our method that should always be taken into account: a 
structure that is very cleanly detected at, say, D=30 kpc 
may go completely undetected if located at D=15 kpc, 
instead. This implies that while significant detections 
of RC peaks in our SCPs are robust and provide accu- 

^^ For each passband, we fitted the error curve derived from 
SDSS photometric errors with exponential functions. Then we 
used the fitted functions to assign the proper average error to each 
synthetic star (according to its magnitude) and we added to each 
synthetic magnitude an error component extracted from a gaussian 
distribution having a equal to the average error. 



rate distances, the significance of non-detections must 
be evaluated with great care, on a case by case basis, 
and, in general cannot be taken as a proof of the ab- 
sence of a given structure (that, for example, may be 
predicted by some model) . One can conceive various dif- 
ferent techniques to mitigate this dependence of sensitiv- 
ity on distance, as, for example, to scale the bin width 
with magnitude to account for the effect described at 
point (1), above, or to transform the magnitudes of any 
color-selected RC candidate into distances, by assuming 
template values for M/ and My, and then to search for 
peaks in Distance Functions instead of SCPs. However 
each of these possible solutions would have an impact 
on the accuracy of the derived distance scale: for this 
reason we prefer to maintain an approach that derives 
differential distances from the direct comparison of truly 
observable quantities of strictly the same nature, i.e. the 
magnitudes of peaks in SCPs that can be determined to 
within a few hundredths of mag, finally providing differ- 
ential distances with accuracies <10%. Other methods to 
trace structures are intrinsically more powerful for other 
purposes (like the detection of structures, for example; 
see Bel06). The technique adopted here is best suited 
for distance measures and we decided to optimize it to 
this task, at least in the present application. 

3. MODELING OBSERVED SCPS AND MEASURING 
DIFFERENTIAL DISTANCES 

In the following we will describe the technique that 
we use to accurately estimate the magnitude (and the 
statistical significance) of peaks in the SCP of on-Stream 
fields. To illustrate it we recur again to the case of the 
Sgr34 field (we will show below also two cases of on- 
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Stream fields). The approach is strictly the same for V 
and / SCPs and it is presented in parallel in the left and 
right columns of panels of Fig. 6: 

1. Upper panels: a possible RC peak is detected in 
both SCPs. The underlying smooth SCP (the 
back/foreground component) is fitted with a sim- 
ple function of the form f{x) = Ae^ + Bx + C, 
once the points in the range enclosing the peak are 
excluded (see also B06a and Correnti et al. 2009). 
This simple procedure is very effective in all the 
cases considered here and it allows a reliable in- 
terpolation of the background in the region of the 
peak. The dotted lines mark the 3,4 and 5cr levels 
above the background, that are computed includ- 
ing both the Poisson noise and the uncertainty in 
the fit of the background component (that is con- 
servatively assumed to be of the same order of the 
Poisson noise, i.e. Vn , where N is the value of the 
model at a given position) . To verify the validity of 
our assumption for running histograms (that have 
non-independent bins), we made several trials (in 
different fields) using "classical" histograms, where 
Poisson statistics clearly apply. We found that the 
statistical significance lines are at the same level in 
both classical and running histograms, confirming 
that our simple assumptions lead to correct results 
(see as another example of the same approach Cor- 
renti et al. 2009). 

2. Middle panels: the model of the background {f(x)) 
is subtracted from the observed SCP. The only sig- 
nificant residual is the RC peak that, in general, 
has a rather symmetric bell shape. The peak is fit- 
ted with a gaussian curve (G(x)) by searching for 
the three parameters of G{x) (mean, a and normal- 
ization factor) that minimize the reduced x^- The 
derived mean is taken as the best estimate for the 
position of the considered peak. As the bins of the 
adopted running histogram are not independent, 
the resulting x^ values can be considered only in 
a relative sense. After several trials on real cases 
we found that increases of x^ by a factor of 2 with 
respect to the best solution (having x^ = Xmin) ^1" 
ways correspond to clearly unsatisfactory fits (see 
Fig. 9, for an example). For this reason we adopt 
the difference between the best-fit mean and the 
mean of the solutions having x^ = "^Xmin ^^ ^ ^O" 
bust estimate of the accuracy of our measures. As a 
sanity check we tested the assumed models against 
the corresponding ordinary histogram (i.e. with 
independent bins) for nine flag=l peaks, comput- 
ing the actual la uncertainty. It turned out that 
our empirically defined error is always within a fac- 
tor from 0.5 to 2 times the statistically correct Icr 
error, thus providing a realistic estimate of the un- 
certainty of our measures. 

3. Lower panels: the global model, obtained by sum- 
ming f{x) and G{x), is compared to the observed 
SCPs. This final form of the overall fit is what we 
will show for all the considered fields in Sect. 4, be- 
low. In Sect. 3.1 we will show that the adopted 
model of the fore/background component of the 



considered SCPs (/(x)) provides an adequate rep- 
resentation of what is observed in actual CFs and 
predicted by current Galactic models. 

The application to the Sgr34 field just described, pro- 
vides also the Zero Points of our differential distance 
scale, i.e. the magnitude of the RC in the main body 
of Sgr, V = 17.82 ± 0.02 and / = 16.87 ± 0.02. As 
a sanity check, we verify if these numbers are compati- 
ble with theoretical stellar model predictions'^"^. Adopt- 
ing the distance modulus (to — M)q = 17.10 ± 0.15 for 
Sgr (Monaco et al. 2004) we obtain M/ = -0.23 ± 0.15 
and My = -1-0.72 ± 0.15. These correspond to ages in 
the range 5-7 Gyr for [Fe/H] = -0.4 and 9-11 Gyr 
for [Fe/H] — —0.7 in the models by Girardi & Salaris 
(2001), in good agreement with all recent estimates of 
the typical age of the bulk of the Sgr stars (see Layden 
& Sarajedini 2000; Monaco et al. 2002, B06a). 

3.1. The SCPs of Control Fields 

To verify empirically that the peaks we interpret as 
due to intersections of the considered los with Stream 
wraps are not due to Galactic structures, we have in- 
spected all the color-selected SCPs of the Control Fields 
described in Sect. 2. The overall conclusion is that there 
is nothing similar to the peaks we observe in the SCP of 
our on-Stream fields in generic Galactic fields at similar 
distances from the plane and the center of the Galaxy. 

In Fig. 10 we show various examples: the SCPs of 
six on-Stream fields (continuous lines) are compared to 
the SCPs of their corresponding CFs (dotted lines, see 
Fig. 1). The best-fit models for the on-Stream SCPs, 
together with the background and the 3(t levels are also 
reported, using the same symbols as in Sect. 3 and Sect. 4 
below. The two SCPs are normalized by the ratio of the 
sampled areas, but any other reasonable normalization 
(for example, by the ratio between the number of stars 
that fall inside our color selection) does not significantly 
change the results. 

The shapes of the SCPs are very similar in the range 
not affected by the peaks associated with the Stream, 
as already observed when we have done the same com- 
parison in the main body (Fig. 6). It is quite clear that 
the strong and well defined peaks observed in on-Stream 
fields are lacking in the SCPs of Control Fields (however 
there is no guarantee that genuine and yet unrecognized 
structures are present also along these los ). It is also re- 
assuring to note that the models for the SCP of the con- 
taminating back/foreground population we have adopted 
for the on-Stream fields provide a good description also 
of the CF SCPs, at least out to F ~ 18.5. Beyond this 
limit it is quite clear that in the on-Stream fields there 
is an additional source of contamination, that has to be 
ascribed to RGB, SGB and MS stars from the Stream 
population itself, as discussed in Sect. 2.1. This provides 
further support to the idea that the adopted approach of 
fitting the back/foreground component directly on on- 
Stream SCPs is the most effective way to get rid of this 

^^ That, however, are quite uncertain and model dependent, in 
an absolute sense, for stars in this evolutionary phase. For instance 
the absolute / magnitude of the peak for a age=6 Gyr, [Fe/H] = 
—0.5 model from the BASTI datasct (shown in Fig. 5), is matched 
by a model of age=1.7 Gyr and [Fe/H] = —0.4, from the set by 
Girardi & Salaris (2001). 
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Fig. 6. — Upper panels: de-reddened running histogram SCPs of color- selected RC candidates for tfie Sgr34 field (continuous line) in V 
and / magnitude (left and right panel, respectively); the dashed line represents the polynomial fit of the background (f(x)). The dotted 
lines mark the 3,4 and 5cr levels above the background, that is computed including both the Poisson noise and the uncertainty in the fit. 
Middle panels: Residuals of the subtraction between the observed SCPs and the fit of the SCPs without the peaks; in red we plotted the 
fit of the peaks, obtained with a gaussian (G{x)). Bottom panels: same as the first ones, with added in red the total fit of the SCPs 
(polynomial for the bkg -|- gaussian for the peaks, f{x) + G{x)). 



kind of self-containation from other species of Stream 
stars, that would not have been possible if we merely 
subtracted the CF SCPs to the on-Stream ones. 

In the lower right panel we present the case of F15A, 
that will be discussed in Sect. 4.1. It is interesting to note 
the close similarity between the two considered SCPs for 
this field, where we do not detect any signal from the 
Stream, and they are therefore expected to be (both) 
dominated by the generic halo/thick disk Galactic pop- 
ulation. 



3.2. Examples of on-Stream fits 

In Fig. 7 and Fig. 8 we show two examples of applica- 
tion to on-Stream fields, the fields F7A and F5A, respec- 
tively. In the first case a broad peak with significance 
above 4a is detected in both the V and / SCPs. The 
derived differential distances with respect to the main 
body are in good agreement, within the uncertainties. 



The f{x) + G{x) model provides an excellent description 
of the observed SCPs. 

Two significant peaks are detected in the SCPs of the 
field F5A (Fig. 8), thus, in this case, we need a model 
with two gaussian components. Both peaks are signifi- 
cantly narrower then that found in field F7A. Neverthe- 
less the model f{x) + Gi{x) + G2{x) provides an excellent 
representation of the observed SCPs. The differential 
distances obtained from the V and / SCPs are in good 
agreement: there is no doubt that we are detecting the 
sam,e structures in both SCPs. 

To place the results shown in Fig. 7 and Fig. 8 into 
the proper context, we plot the positions of the detected 
peaks into the XQ^SgrY^^sgr plane of the heliocentric 
Sgr coordinates as defined by Majewski et al. (2003), 
in Fig. 11. This plane is defined to coincide with the 
plane of the orbit of Sgr, hence the Stream is expected 
to be confined within a few kpc about it. We compare 
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Fig. 7. — Upper panels: de-reddened running histogram SCPs of 
color-selected RC candidates for a Stream field (field F7A, continu- 
ous line) in V and / magnitude (left and right panel, respectively); 
the dashed line represents the polynomial fit of the background 
(/(x)). The dotted lines mark the 3,4 and Scr levels above the 
background, that is computed including both the Poisson noise 
and the uncertainty in the fit. Middle panels: Residuals of the 
subtraction between the observed SCPs and the fit of the SCPs 
without the peaks; in red we plotted the fit of the peaks, obtained 
with a gaussian {G{x)). Bottom panels: same as the first ones, 
with added in red the total fit of the SCPs (polynomial for the bkg 
-I- gaussian for the peaks, f{x) + G{x)). 




Fig. 8.— As Fig. 7, but for a difll'erent field (field F5A). The 
method is the same although this field shows more than one peak. 

our detections with one of the N-body models of the tidal 
disruption of Sgr by Law et al. (2005). In particular we 
plot in Fig. 11 the results of the evolution of the N-body 
model of Sgr within a Galactic DM halo of prolate shape 
(flatness q=1.25, see Law et al. 2005; Johnston et al. 
2005, for further details on the models). To compare ob- 
servations and model in a consistent way we transformed 
our relative distances into absolute values by adopting 
the same distance modulus for Sgr as Law et al. (2005), 
i.e. (m-M)o = 16.9 (Mateo et al. 1995a). The points of 
the model that are encountered by the considered F.o.V.s 
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Fig. 9. — Upper panels: de-reddened running histogram SCPs, 
zoomed in the region of the peak(s) , of color-selected RC candidates 
for the Stream field F5A (left panel) and F7A (right panel). The 
meaning of the lines is the same as Fig. 7, the red line represents 
the global model, f{x) + G{x), with the best fit value of G{x) mean 
(value for which x^ = Xmin)- Middle and lower panels: as upper 
ones, with the exception that the values of the G{x) means are 
those that have x^ = ^Xmin i X^ = —'^Xmiw middle panels and 

^ — +'^Xmin' lower panels, respectively). It is clearly visible that 




X 

in these last two cases the fit is totally unsatisfactory. 




Fig. 10. — Comparison between a sample of the on-Stream SCPs 
studied in this paper (continuous histogram; best-fit model in red) 
and the SCP of the corresponding Control Fields. 

along the los (F7A and F5A, from left to right, respec- 
tively) are plotted as heavier dots. 

Taking the considered model as a realistic representa- 
tion of the actual Sgr relic (a very reasonable assumption, 
in first approximation; Law et al. 2005), it is clear that 
any los around the considered plane would cross one or 
more different wraps of the Stream, at different distances 
(see Fig. 11). The peak from F7A and the most distant 
peak from F5A seem to match a distant portion of the 
leading arm. The nearest peak from F5A matches very 
well with a wrap of the trailing arm that appears narrow 
and well defined and that is crossed nearly perpendicu- 
lar by the considered los . According to the considered 
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Fig. 11. — N-body model of the tidal disruption of Sgr by Law 
et al. (2005, prolate Dark Matter halo case) in heliocentric Sgr 
coordinates (see Majewski et al. 2003). The heavily printed dots 
are the subset of the model particles that are enclosed in the two 
observed fields considered here (F5A and F7A). The observed po- 
sitions of the Stream in these fields, as estimated from the position 
of the RC peaks in our SCPs, are plotted as filled (red) circles. 

model, both los should also cross a nearby wrap at a 
distance not enclosed in our range of sensitivity, that is 
delimited by the two dashed circles in Fig. 11. No wrap 
is expected to lie outside D=60 kpc in the region sampled 
by our fields. 

Both the more distant F5A detection and the single 
F7A detection occur in regions where the model predicts 
the confluence and crossing of different wraps. At a first 
glance to Fig. 11, it may appear that the constraining 
power of a single "mean position" of a Stream wrap, as 
derived with our method, is not sufficient to describe the 
complex structure of the Sgr remnant along a given los . 
In Fig. 12 we compare the observed RC peaks of F5A and 
F7A (and their best-fit gaussian models) with the peaks 
derived from the N-body model shown in Fig. 11 and 
from the oblate-halo model from the same set (Law et al. 
2005) by (a) selecting the model particles encountered by 
the FoV cone, (b) assigning to each of them the absolute 
magnitude of the RC {M^^ = -^0.72 and Mf^ = -0.23, 
according to Sect. 3) and deriving their apparent mag- 
nitude according to their distance, (c) adding gaussian 
photometric errors as a function of the apparent magni- 
tude similar to the observed ones, and (d) producing the 
running histogram of the derived magnitudes with the 
same settings adopted for the observed SCPs. To make 
easier the comparisons shown in Fig. 12, the synthetic 
SCPs have been multiplied by an arbitrary normalization 
factor. The qualitative resemblance of the observed and 
predicted structures for the prolate-halo model is strik- 
ing. On the other hand the oblate-halo model is clearly 
unable to reproduce the observations, even in terms of 
number of Stream wraps encountered by the considered 
los . 

A more thorough comparison of our observations with 
theoretical models of the disruption of Sgr will be pre- 
sented in Sect. 6. Here we are just interested to demon- 
strate that our method allows a detailed comparison be- 
tween models and observations not only in terms of mean 




Fig. 12. — Comparison between the observed field-subtracted 
RC SCPs from F5A (upper panels) and F7A (lower panels) and the 
predictions of the N-body models by Law et al. (2005) in the same 
directions (left panels: prolate-halo model, the same as Fig. 11; 
right panels: oblate halo model). The observed SCPs are the con- 
tinuous running histograms, the red continuous curves are the best- 
fit gaussians, the dotted running histograms are the predictions of 
the models, including the photometric errors. The histograms from 
the N-body models have been multiplied by arbitrary normaliza- 
tion factors (2x and 3x, for the left and right panels, respectively) 
to make the plots more readable. 

distances, but also in terms of the actual shape of the 
structures along the los (see Fig. 12, left panels, in partic- 
ular) . In other words a fully successful model of the Sgr 
Stream must reproduce the correct position and shape of 
the observed peaks: this provides the opportunity for a 
fruitful detailed comparison between models and obser- 
vations also in regions were different Stream wraps cross 
each other. In the present study we provide the position 
and the FWHM of the peaks (as the most basic shape 
parameter, see Sect. 3.5) but anyone interested in more 
detailed comparisons can easily reproduce our results and 
obtain plots like Fig. 12. 

In Sect. 4 the SCPs of all the considered on-Stream 
fields, with the associated detections will be presented 
and briefiy discussed. In agreement with the qualitative 
predictions of the model shown in Fig. 11, in most cases 
we will detect two peaks at different distances. 

3.3. Classification of the detections 

We divided the detections in three categories, assigning 
a flag to each of them, according to the following criteria: 

• flag ~ 1: peaks that are above the 3cr thresh- 
old both in the V and / SCPs, and having the 
same A(Mag) in both passbands, within the er- 
rors. These are called primary peaks. 

• flag — 2: peaks having the same A{Mag) in both 
passbands, within the errors, but reaching the 3a 
threshold only in one of the two SCPs. In all 
the flag = 2 cases described in the following, the 
weaker peak is always just below 3cr. These are 
called secondary peaks. 

• flag — 3: clearly visible peaks having the same 
A{Mag) in both passbands, within the errors, but 
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not reaching the 3(7 threshold in both the SCPs. 
These correspond to uncertain detections that we 
report just for completeness. In some case a weak 
peak in the SCP in one band is tentatively identi- 
fied as it corresponds to a stronger peak in the other 
SCP. These peaks are called also tertiary peaks. 

The observed SCPs and the detected peaks will be 
briefly described and discussed in Sect. 4. 

3.4. Number of RC stars associated to each peak 

Our modeling of the observed SCP automatically pro- 
vides also an estimate of the total number of stars asso- 
ciated to any given peak. This gives a useful additional 
constraint for theoretical models; we will illustrate this 
possibility with an example in Sect. 6.2. In Tab. 2 we 
provide the number of stars associated with a given peak 
normalized to an area of 25 deg^ (quite similar to the 
actual area of our fields). This number is the weighted 
mean of the estimates obtained from the V and / SCPs, 
where the assumed error on the estimate from each SCP 
is just the square root of the observed number (hence it 
should be considered as a lower limit to the real error). 

Adopting the Sgr34 field as a baseline (having 
'. 1500 stars/deg^, and muy ~ 25.5 
mag/arcsec^; Monaco et al. 2005b), and translating our 
on-Stream estimates into stars/deg^ units {N^'~' , for the 
field i) we can transform our numbers into surface bright- 
ness, according to the same formula used in Bellazzini 
et al. (2006b): 



prRC 



MV,J = My,Sgr34 - 2.5log{ 



N. 



RC 



atRC 
^^Sgr34 



) + A(m-M)o (4) 



derived from Renzini (1998), where A(to — M)q = (m 



MY, - (m - Mr, 



Sgr3i 



AV — Al. Using the equation 



4 we find that the V surface brightness of the portions 
of the Sgr Stream studied in the present paper range 
between 30.6 mag/arcsec^ and 33.6 mag/arcsec^, quite 
typical of tidal tails and in good agreement with previ- 
ous results (see, for example, Bellazzini et al. 2003b, and 
references therein). 

It is important to recall that the measured densities 
refer only to RC stars: in the presence of a population 
gradient (as is likely the case in the Stream, Bellazzini 
et al. 2006c; Chou et al. 2007) they would trace differ- 
ent fractions of the total stellar content at different po- 
sitions along the Stream. Analogously, the derived sur- 
face brightness estimates have been rescaled assuming 
the stellar mix of Sgr34 for all the considered portions 
of the Stream. For this reason these estimates should be 
considered with caution: given the sense of the gradient 
it is expected that they provide lower (faint) limits when 
converted into luminosity or surface brightness. In this 
context, it is interesting to note that if we convert our 
surface brightness into the same density units (L0/kpc) 
adopted by Niedcrste-Ostholt et al. (2010), we find that 
our results are fully compatible with the trend of den- 
sity as a function of RA derived by these authors (see 
their Fig. 7), for both branches. On the other hand, our 
densities are lower than theirs by a factor of ^ 4 — 5. It 
is reasonable to assume that part of the difference may 
be accounted for by the effect of the population gradient 
described above. 



3.5. Depth along the los of the Stream wraps 

Fig. 12 shows that the observed RC peaks contain also 
valuable information on the characteristic size of the sec- 
tion of the Stream branches crossed by our los, as peaks 
at similar distances display different widths. To obtain a 
quantitative estimate of the linear width along the los of 
the structures identified here, we recurred to the syn- 
thetic population described in Sect. 2.3. In particular 
we tried to reproduce the models of the observed peaks^** 
with smoothed histograms of the synthetic RC popula- 
tion, properly including the effects of photometric errors. 
As done in Sect. 2.3, we assign a distance along the los to 
each star of the synthetic RC population according to a 
gaussian distribution having the same mean and normal- 
ization as the observed peak, and we search for the value 
of ad giving the best match between the two models. In 
Tab. 2 we report the Full Width at Half Maximum (in 
kpc) of the adopted distributions, FWIIM=2.35crd- The 
best match is found by minimizing x^ i typical uncertain- 
ties are ^ 20%. The adopted procedure gets rid of the 
effects of the distance on the width of the SCP peaks 
discussed in Sect. 2.3. 

Since the synthetic population that we adopt is strictly 
single-age and single-metallicity, the intrinsic luminos- 
ity width of its RC should be smaller than the actual 
width of the RC of Sgr, as the latter hosts stars span- 
ning a range of ages and metallicities (B06a, Girardi & 
Salaris 2001; Siegel et al. 2007; Bellazzini et al. 2008). 
For this reason the FWHM values we obtain in this way 
must be considered as strong upper limits to the real 
values. Moreover, it has to be recalled that we report 
FWHM along a given los , that may have various in- 
cidence angles with respect to the encountered Stream 
wraps. Applying the method to the SgrS4 field we ob- 
tain FWHM ~ 3 kpc, not too far from the minor-axis 
FWHM in the plane of the sky as obtained from the 
best-fit King (1962) model by Majewski et al. (2003), 
i.e. FWHM ~ 1.1 kpc, in particular if we take into 
account that the los toward the core of Sgr is (likely) 
not exactly perpendicular to the major axis of the dwarf 
galaxy (see Fig. 11). Based on this test, it is reasonable 
to assume that our FWHM overestimates the true values 
by a factor of > 2. 

In any case, the ratio between the FWHM of two dif- 
ferent los /locations in the Stream, or, equivalently, the 
differential trend of the FWHM as a function of orbital 
azimuth along a given Stream wrap, can be directly com- 
pared to the predictions of theoretical models of the dis- 
ruption of Sgr. 

4. ON-STREAM DETECTIONS 

In this section we present all the SCPs obtained from 
each analyzed field of branches A and B; we plot the 
SCPs in the V- and /-bands in Fig. 15 for branch A, and 
Fig. 16 for branch B. Together with the observed SCPs, 
we plot also the background model (/(x), dashed lines), 
the threshold limits for the detections (dotted lines, re- 
spectively 3,4, 5cr), which as mentioned previously, in- 
clude both the Poisson noise and the uncertainty in the 

^"^ In this way we avoid any problem associated with the (partial) 
overlap of observed peaks; adjacent overlapping peaks are disentan- 
gled by our models as Point Spread Function -fitting photometric 
packages disentangle the fluxes from two partially overlapping stars 
on an image. 
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fit, and the global model that fits the observed SCPs 
(/(x) + G(x), red lines). Each field is labeled according 
to the names assigned in Table 1 and with its Galactic 
coordinates (l,b). 

In summary, we detect 26 primary (flag=l) peaks, 10 
secondary (flag=2), and 14 tertiary (flag=3) peaks. Most 
of the considered SCPs show two significant peaks, cor- 
responding to subsequent crossings of different wraps of 
the Stream along the los . The trend of peak distance 
as a function of Sgr longitude (Asgr) shown in Fig. 13 
can be useful to better understand the morphology of 
the various SCPs presented below. Most primary peaks 
appear to trace a wrap of the leading arm whose dis- 
tance from the Sun steadily decreases from D ~ 45 kpc 
at Asgr ^ 290° to £» ~ 20 kpc at Asgr ^ 230°. Both 
primary and secondary peaks trace a more nearby fila- 
mentary structure at constant distance D ~ 25 kpc, from 
Asgr — 290° to Asgr — 260°, that then bends toward 
larger distances, reaching D ~ 40 kpc at Asgr — 230°. 
This feature is tentatively identified as a wrap of the 
trailing arm (see Sect. 6). The two wraps cross at 
Asgr ^ 245° (see Sect. 6). Some secondary and ter- 
tiary peaks seem to trace more feeble distant or nearby 
wraps (see Sect. 6 for a deeper discussion). The compar- 
ison with the considered model suggests that most of the 
detected peaks can be associated with the Sgr Stream. 
The tertiary peak at Asgr ^ 263° and _D ~ 18 kpc, and 
the primary peak at Asgr ~ 280° and D ~ 19.5 kpc, are 
possibly associated to other overdensities in the Virgo 
constellation (see Juric et al. 2008; Duffau et al. 2006; 
Newberg et al. 2007, and references therein), as discussed 
in some detail in Sect. 6. 

Before proceeding in the description of the various de- 
tections, we anticipate that the differential distance mod- 
uli (Ay = VJiQ^'^ — Vj^^ ; the analogous definition be- 
ing valid also for A/), reported in Tab. 2, obtained from 
primary and secondary peaks detected in the V and / 
SCPs are in excellent agreement, as shown in Fig. 14. 
This confirms the reliability of our detections and dis- 
tance estimates. For this reason, from Sect. 5 on and in 
Fig. 13 we adopt the mean of /S.V and A/ as our final 
differential distance moduli estimates. 

4.1. Branch A detections 

In branch A we analyzed 15 fields, the correspond- 
ing observed SCPs and the adopted best-fit models are 
shown in Fig. 15. We obtained a total of 24 peak detec- 
tions, with the following classification: 10 primary peaks, 
5 secondary peaks and 9 tertiary peaks. The SCPs of 
the first five fields (from FIA to F5A, upper left panel 
and first two rows of the upper right panel of Fig. 15) 
display a common general behaviour: they present two 
main peaks, the one at fainter magnitudes always be- 
ing the strongest (a primary peak in all cases), while 
the brighter ones are wider and span all the classes from 
flag=l to fiag=3, depending on the specific field. It is 
quite clear that this series of peaks traces two coherent 
structures placed at different distances along the los . 
The /-band SCP of F4A may suggest a splitting of the 
brighter/weaker peak into two separate components: we 
consider this interpretation as unlikely, nevertheless the 
result obtained with a three peaks model is briefly dis- 
cussed in Sect. 4.3. The only exception is a primary peak 
detected ai V ^ 16.5 in FIA: this likely corresponds to 



the nearest wrap of the Stream that emerges from the 
d < 12 kpc circle (where our method is blind), which is 
the most Eastern of all the considered los (see Sect. 6). 

The SCPs of F6A (second upper right panel of Fig. 15) 
present an overall structure similar to those described 
above. However we identified additional (fainter) peaks, 
and finally we adopted a three peak solution, whose valid- 
ity is confirmed by the inspection of SCPs obtained with 
reduced bin width (i.e., higher resolution^^; see Sect. 4.3 
for an alternative). The newly- resolved third peak, at 
magnitude F ~ 19 (/ ^ 18), corresponds to a distance 
I? ~ 42 kpc at Asgr — 265°; this detection seems related 
to a very distant wrap of the leading arm (see for exam- 
ple Fig. 22). F6A is the only branch-A field in which we 
detect a peak related to this distant wrap of the leading 
arm, that was observed also by Bel06 (see Sect. 5); the 
same structure is detected in branch B along the same 
los , as well as along an additional one (F8B). 

F7A is one of two fields in branch A that presents 
only one detection: SCPs (third upper right panels of 
Fig. 15) show a single, very prominent primary peak. As 
discussed above, this los intercepts a region where two 
or three wraps of the Stream cross each other. Fig. 12 
shows that the presence of a single peak is nevertheless 
consistent with model predictions. In both the SCPs of 
the adjacent field F8A, a remarkably weak peak appears 
at similar position as in F7A, hence we obtain only a ten- 
tative flag=3 detection. We have no convincing explana- 
tion for the weakness of the peak detected in this field: 
it may be related to the complex structure of the vari- 
ous Stream wraps or to a local dip in the density along 
the Stream. However, the derived distance is in good 
agreement with the trend observed in others Branch-A 
fields. In the SCPs of F9A field we identify again two 
peaks, the faintest one being very prominent and wide; 
also this los intercepts a region of crossing wraps, thus 
superposed structures may contribute to the production 
of a strong and remarkably wide primary peak. The 
weaker/brighter peak is more interesting: it is clearly 
identified in both SCPs, even if below the 3ct threshold, 
at y ~ 17.20 (/ ^ 16.30): as discussed later in Sect. 6 
this feature may trace a near wrap of the Stream that 
was never detected before. 

SCPs of fields from FlOA to F12A show two peaks 
at similar positions, with a remarkable variety of abso- 
lute and relative strengths. This may refiect the highly 
structured morphology that is suggested by models in 
this region (see Sect. 6). FllA is crossed by the Orphan 
Stream (Bel06, Belokurov et al. 2007). While the dis- 
tance of this structure (~ 30 kpc toward this los ) does 
not match with the detected peaks we cannot rule out 
some contamination from Orphan Stream stars in this 
field. 

We did not find any convincing signal in F13A, F14A 
and F15A; the overall shape of the SCPs appear quite 
different from the other cases and, in the case of F15A, 
the polynomial model did not provide a satisfactory fit 
to the background population. In particular, the SCPs 
present a strong excess at bright magnitudes {V < 18, 
/ < 17) with respect to those in the previously-discussed 
fields, such that they appear flat or even decreasing with 

^^ We note that this is the only case in which a change in the 
bin width produced a change in the interpretation of the SCPs. 
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Fig. 13. — Distribution of the primary (red filled circles - branch A -, and blue filled squares - branch B), secondary (red open circles 
- branch A -, and blue open squares - branch B) and tertiary detections (red starred symbols - branch A, and blue stars - branch B) in 
the Aggj. vs. heliocentric distance plane (a true distance modulus of 16.90 has been adopted here). The horizontal dashed lines enclose 
the range of sensitivity of our method. The prolate-halo N-body model by Law et al. (2005) is also reported (small dots) as an aid for the 
interpretation of the plot. The heavier dots are those enclosed in the cones of the considered FoVs. 



increasing magnitude. These are the fields at the low- 
est galactic latitude, hence we attribute these features 
to contamination by (relatively) nearby stars from the 
Galactic thick disk (and, possibly, the Monoceros struc- 
ture, see Fig. 1 of Bel06) that overwhelms the signal from 
the Sgr Stream RC. This seems confirmed by the com- 
parison with the corresponding CFs (the case of F15A 
is shown in Fig. 10, in Sect. 3.1), that display SCPs es- 
sentially indistinguishable from those of the on-Stream 
fields. This implies that the adopted technique can be 
used successfully only at large distances from the Galac- 
tic plane. Given the above reasons we preferred not 
to consider for further analysis the possible peaks at 
V(I)- 16(15) and V(I)~ 17.5(16.5) in F14A. 

4.2. Branch B detections 



In branch B we analyzed 13 fields, obtaining a total 
number of 26 detections, with the following classifica- 
tion: 16 primary, 5 secondary and 5 tertiary peaks. The 
main structures found in branch A are mirrored also in 
branch B, as clearly shown in Fig. 13. In all the fields 
(except F7B) we detect at least two peaks; in two cases 
(F8B and F12B) we also detect a third peak; in FSB this 
is likely tracing the more distant Stream wrap running 
nearly parallel to the main wrap of the leading arm (see 
above and Fig. 13); in another case (F5B), in addition 
to three peaks analogous to those in FSB, we found an 
additional nearby peak (see Sect. 6 for a discussion). An 
alternative interpretation for the SCPs of F12B is pre- 
sented in Sect. 4.3. Quite suprisingly, the SCPs of F6B 
appear completely smooth and featureless. In this case 
we were not able to find an explanation for this behavior 
(but see above for the discussion of the similar case of 
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TABLE 2 
Distance, FWHM and density for the detected RC peaks. 
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F8A). 

In analogy with F7A (and F8A) F7B is the only case 
of branch B SCPs fitted with a single peak model. The 
morphology of the /-band SCP and the comparison with 
the adjacent FSB field suggest that two, or possibly three 
peaks may be merged together in this SCP. However we 
were unable to resolve the peak into separate components 
even in SCPs with smaller bin width (as for the case of 
F6A). We caution the reader that this primary detection 
is likely concealing significant — but as yet undetected 
— substructure. 

4.3. A few special cases 

There are a few cases in which the observed SCPs do 
not provide unequivocal indications for the model to be 
adopted, in particular concerning the number of Gi{x) 
functions to be included in the model, i.e. the number 
of detected peaks. F7B, briefiy discussed above, is the 
only case in which we feel that the observed peak is due to 
the merging of two (or, more likely, three) adjacent peaks 
that we cannot resolve. In Fig. 17 we present acceptable 



alternative models (with respect to the solutions shown 
in Figs. 15 and 16, and listed in Tab. 2) for the three cases 
in which our preference for the adopted models (Tab. 2) 
is only marginal, and is also supported by the continuity 
within a large scale structure (a Stream wrap). In Tab. 3 
we report the corresponding alternative solutions, that 
can be replaced with those of Tab. 2 by those readers 
who may use our values to constrain models of Sgr, if 
they judge them more appropriate, for some reason. 

4.4. Intra-Branch fields 

In the present analysis we do not consider the structure 
of the Stream in the Dec direction. We fully adopted the 
view of Belokurov et al. (2006), where the leading arm 
of the Stream as seen from TO stars in the SDSS bifur- 
cates into branch A and branch B around RA=220° and 
the separation between the branches increases with de- 
creasing RA. We proceeded to a basic verification of this 
scenario by looking at the SCPs of a few Intra-Branch 
(I) fields (not shown here, for brevity), located at inter- 
mediate Dec with respect to the A and B fields F5, F7, 
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Fig. 14. — Comparison between the differential distance moduli 
obtained from peaks in the V and I SCPs for primary (upper panel) 
and secondary (lower panel) peaks. 

FIO and F12. In agreement with the results of Belokurov 
et al. (2006), we find that the SCPs of F5I and F7I mimic 
the structure of the SCPs of the corresponding A and B 
fields, showing peaks at the same position and with sim- 
ilar shape, but weaker than in the on-Stream fields (i.e. 
tracing a lower stellar density). In the SCPs of FIOI the 
peaks seen in the A and B SCPs are just barely visible 
and they completely disappear in F12I. Hence, these lim- 
ited set of tests confirm the reality and the morphology 
of the Stream bifurcation as observed by Belokurov et al. 
(2006). 

4.5. The color of the RC peaks 

The color of the RC peak is known to be very sensitive 
to metallicity and weakly sensitive to age, in the range of 
ages relevant for the present study (4-12 Gyr, Girardi & 
Salaris 2001). As our procedure of independent peak de- 
tections in V- and /-band SCPs automatically provides 
the colors of the RC peaks, it is worth checking if there 
is any hint of a color (metallicity) gradient along the 
Stream. In Fig. 18 the colors of the observed peaks are 
compared to the theoretical models by Girardi & Salaris 
(2001). All of the peak detections shown in Fig. 18 are 
compatible with having the same color within the un- 
certainties (that are quite large for some non-primary 
peaks). It is interesting to note that the large majority of 
points cluster around the [M/ H] = —0.7 model, in good 
agreement with the results by Bellazzini et al. (2006a) 



and Carrell & Wilhelm (2010) (see also Sect. 2.1). No 
significant trend of color (metallicity) of the RC popula- 
tion with orbital azimuth is apparent and the few points 
showing the larger color difference from the mean locus 
are always among those having the most uncertain color 
estimates. 

5. COMPARISON WITH PREVIOUS ANALYSES 

Before discussing in detail the comparison between 
our distance estimates and the findings from previous 
works, it is worth considering the difference between the 
performances of the various adopted tracers. The in- 
trinsic stability (and ubiquity along the Stream) of our 
standard candle (RC stars), the adopted analysis, best 
suited for the detection and location of RC peaks, and 
the purely differential nature (Stream vs. main body) 
of our measures, make our distance estimates the most 
comprehensive, accurate and homogeneous set publicly 
available (even if limited to the region of sensitivity de- 
scribed above). The uncertainties associated with our 
estimates are lower than any previous work, with a typ- 
ical values of < 5%, raising to < 10% in the worst cases. 
For example, Majewski et al. (2003) reports that the 
characteristic uncertainty of their photometric parallaxes 
based on M-giant is ~ 20%; Martin et al. (2004) showed 
that uncertainties in the age/metallicity of the considered 
populations may lead to systematics of order r^ 30% in 
the distance scale based on M giants. F stars (assumed 
to be TO stars of Sgr) proved to be an excellent mean 
to trace even very feeble substructures (Belokurov et al. 
2006; Newberg et al. 2007). However the assumption of a 
common absolute magnitude for all color-selected F stars 
implies large uncertainties, as these stars span a range of 
luminosities much larger than RC stars. For example, if 
we consider the distribution in V magnitude of (a) the 
RC selected with our color window, and (b) the MSTO 
stars selected in color as done by Bel06 (and limited to 
V > 20.0) in the photometry of the Sgr34 field, we found 
two obvious single peaked distributions, but while the 
FWHM of the RC peak is ~ 0.3 mag, the MSTO star 
peak has FWHM ~ 2 mag. Indeed, Cole et al. (2008), 
in their pilot project on stripe 82, showed that the as- 
sumption of a fixed magnitude for these stars may lead to 
very large errors. Blue Horizontal Branch stars are eas- 
ier to select against the Galactic fore/background, but 
are rarer than RC stars. Moreover, even if selected in a 
color range where the Horizontal Branch is really nearly 
horizontal, the distribution in magnitude of these stars 
is not expected to be as clearly peaked as the RC (see 
B06c). In this sense, the Sub Giant Branch (SGB), used 
by Bel06 and Keller (2009b), is more promising, as it is a 
very narrow feature in CMDs of metal rich populations. 
However it should be much more sensitive to metallicity 
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Fig. 15. — Fits of the observed SCPs (in V and /) for fields on Branch A of the Sgr Stream. The numbers in parentheses are the Galactic 
longitude and latitude of the center of the field, in degrees. The meaning of the symbols is the same as in Figg. 7 and 8, above. No global 
fit has been attempted for SCPs that did not show significant RC peaks (F13A, F14A, F15A). 



and age variations than the RC (see, for example, B06a, 
and references and discussion therein), and being much 
(intrinsicaUy) fainter, its use is hmited to a lower dis- 
tance range, for any given dataset. Finally, RR Lyrae 
stars (Ivezic et al. 2000a; Vivas & Zinn 2006; Prior et al. 
2009a; Keller et al. 2008) can provide distances with even 
superior accuracy with respect to our method; well sam- 
pled light curves can also give indications on physical 
properties of individual stars (metallicity, for example) 
that cannot be obtained from RC stars. However RR 
Lyrae are (likely) less frequent than RC stars over most 
of the Stream extension and, above all, they need time 
series information to be safely identified and to obtain a 
reliable apparent magnitude averaged over the pulsation 
period: for this reason the available data cover a much 
smaller region of the sky with respect to generic "single 
epoch" standard candles. 

5.1. Comparison with specific detections 



Yanny et al. (2000) were the first to interpret a stel- 
lar over-density in the halo as possibly due to the Sgr 
Stream. In the first available (equatorial) stripe of the 
SDSS they identified an excess of A-type stars around 
A© ~ 295°, adjacent to our field FIA. The heliocentric 
distance inferred is of 48 kpc, in good agreement with our 
estimate for the main wrap of the leading arm in this di- 
rection {D ~ 45 kpc at Asgr — 290°). This result was 
later confirmed by the more thorough study by Newberg 
et al. (2002), that used F stars as main tracers. Although 
they do not comment on it, the Yanny et al. (2000) data 
also showed an excess of A-type stars less than 20 kpc 
away along the same las (see their Fig. 18 and 19). This 
may be more easily identified with the constant-distance 
coherent structure we see at d ~ 25 kpc than to the 
nearest wrap that we (possibly) detect at A ~ 287° and 
d ~ 13 kpc. 

A similar detection of two density enhancements to- 
ward the Northern Loop was reported by Ivezic et al. 
(2000a), from the study of RR Lyrae in the same SDSS 
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Fig. 16. — The same as Fig. 15 but for fields on Branch B of the Sgr Stream. No global fit has been attempted for SCPs that did not 
show significant RC peaks (F6B, FllB). 

stripe studied by Yanny et al. (2000), and by Vivas & 
Zinn (2006), also using RR Lyrae from the QUEST RR 
survey, which explored nearly the same region of sky 
(Aq ~ 270° — 290°). Both studies comment primarily on 
an excess of RR Lyrae stars at 45-50 kpc (corresponding 
to the main wrap of the leading arm); however a struc- 
ture around '^ 20 kpc is also noted. 

Majewski et al. (2003), provided a clear panoramic 
view of the Sgr Stream using M giants as standard can- 
dles; they were able to trace very neatly the trailing tail 
all over the Southern Galactic hemisphere, as well as part 
of the leading arm closer to the main body of the galaxy, 
up to RA ~ 190°. They report two cases of M giants ex- 
cess along the los in common with the present analysis. 
The most evident at a distance D ~ 45 kpc, compatible 
with our estimates, and the other one, less pronounced, 
at a distance D '^ 20 — 25 kpc, for which the interpre- 
tation is not so clear as in the case of A stars and RR 
Lyrae detections. 

All the detections mentioned above, as well as others 




Fig. 17.— Fits of the observed SCPs (in V and /) for fields of 
the branch A and B of Sgr Stream for which we present acceptable 
alternative models, with respect to the solutions shown in Figs. 15 
and Figs. 16. The numbers in parentheses are the Galactic longi- 
tude and latitude of the center of the field, in degrees. The meaning 
of the symbols is the same as in Figg. 7 and 8, above. 
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Fig. 18. — De-reddened color of the observed RC peaks as a 
function of A^Sgr- The observed points are compared with the 
predictions of the Girardi & Salaris (2001) theoretical models for 
different metallicities, plotted as a function of age of the population 
in the range 4-12 Gyr. The scale on the lower horizontal axis (age) 
refers to models, that on the upper horizontal axis {A^sgr) refers to 
observed points. 

toward specific directions, also not included in the range 
considered here (Martinez-Delgado et al. 2001, 2004; Bel- 
lazzini et al. 2003b; Vivas & Zinn 2006), are collected and 
reported in Fig. 17 of M03. This figure, as well as Fig. 19 
in Law & Majewski (2010), clearly illustrates how it may 
be difficult and misleading to put results from different 
sources (and on different distance scales) all together. 
In this sense, it is more fruitful to compare our results 
with other data sets providing homogeneous distance es- 
timates for significant portions of some wrap in common 
with those detected here. 

For instance, Bel06, who used A-F dwarf stars from 
the SDSS to trace the Stream, detected a distant gra- 
dient along the main wrap of the leading arm that is 
in good agreement with our results (for both branches). 
More interestingly, Bel06 found a double detection in a 
few branch A los (from F5A to F7A): in addition to the 
main wrap of the leading arm, they found also a more 
distant structure, ~ 15 kpc behind. This finding is also in 
excellent agreement with our results (see Sect. 4.1). The 
only difference is that we detect this structure, at similar 
distance, also in the corresponding branch B fields. 

Newberg et al. (2007) investigated the relationship be- 
tween several previously-identified substructures in the 
direction of Virgo and the Sgr Stream using imaging and 
spectroscopic observations of F stars and BHB stars from 
SDSS and SEGUE. In their Tab. 1, they reported the de- 
tections associated to the Sgr Stream, providing also es- 
timates of the distance of these structures. This allowed 
us to perform the direct comparison with our results that 
is presented in Fig. 19. The agreement for the structure 
detected in both studies (main wrap of the leading arm) 
is very good, both for branch A and B. 

Nicderstc-Ostholt et al. (2010) investigated the lead- 
ing arm of Sgr Stream in the same region of the "field 
of stream" analyzed by Bel06, using BHB candidates 
from the SDSS. The distances to BHBs are calculated 
assuming an absolute magnitude Mg — 0.7 and are in 
good agreement with the results of Bel06. In Fig. 19 we 
compare our estimates with those obtained by Niederste- 
Ostholt et al. (2010), reported in their Tab. 1, in the 
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Fig. 19. — Comparison between our distance estimates (same 
symbols as in Fig. 13) and those by Newberg et al. (2007, open 
triangles) and by Niederste-Ostholt et al. (2010, open diamonds). 
The match between the different indicators is remarkably good, 
for the portions of the leading arm considered by the three studies. 
The trailing arm is detected only in the present study. 

region of sky going from RA > 160° to RA < 220°. 
Our distances are calculated using a true distance mod- 
ulus of (m — M)o = 17.10, roughly the same adopted by 
Niederste-Ostholt et al. (2010) (the g magnitude of the 
BHB in the main body is g ~ 17.80). In all the regions 
where the data-sets overlap the match is very good, both 
for branch A and B detections. The trailing arm of Sgr 
Stream is not sampled by Niederste-Ostholt et al. (2010) 
that concentrated their analysis on the detection of the 
leading arm and on an accurate distance estimate for this 
wrap of the Stream. 

In conclusion, the overall agreement with previous de- 
tections of the leading arm is very good. The situation for 
the other coherent structures detected here is more dif- 
ficult to judge; in our view the only firm conclusion that 
can be drawn is that several independent studies found 
evidence for some structures located in front of the main 
wrap of the leading arm, in the considered range of A. It 
is unclear if some of these detections can be associated 
with the constant distant (putative) wrap of the trail- 
ing arm detected here or to even more nearby wraps. In 
this sense it is interesting to note that a similar coherent 
structure, at a similar distance, is detected also by Keller 
(2009b), using SGB stars (see his Fig. 7). 

6. COMPARISON WITH MODELS 

As soon as it was realized that Sgr was likely undergo- 
ing tidal disruption, several authors attempted to model 
the process by means of N-body simulations, to establish 
the plausibility of proposed models and to infer the prop- 
erties (mass, orbit) of the original system (Velazquez & 
White 1995; Johnston et al. 1995; Ibata et al. 1997; Edel- 
sohn & Elmegreen 1997; Ibata & Lewis 1998; Gomez- 
Flechoso et al. 1999; Johnston et al. 1999; Jiang & Bin- 
ncy 2000; Helmi & White 2001). It is interesting to note 
that Velazquez & White (1995) were able to provide es- 
timates of perigalactic and apogalactic distances and or- 
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bital period remarkably similar to those obtained in the 
most recent studies, just one year after the discovery of 
Sgr {Rperi — 10 kpc, Rapo — 52 kpc and Port — 0.76 
Gyr, to compare, for instance, with Rperi — 15 kpc, 
Rapo — 60 kpc and Port — 0.85 Gyr, from Law et al. 
2005). The possible role of Sgr in the formation of the 
Galactic Disk warp was studied by Ibata & Razoumov 
(1998) and Bailin (2003). 

However, since the Sgr Stream appears as a 
remarkably-coherent structure crossing a large part of 
the Galactic halo on a nearly-polar orbit, it seems the 
ideal tracer to study the overall shape and the degree 
of clumping of the Galactic halo as a whole. For this 
reason, the most recent N-body modeling efforts have fo- 
cused on constraining the shape of the DM halo of the 
Milky Way (Ibata et al. 2001b; Martinez-Delgado et al. 
2004; Helmi 2004; Law et al. 2005; Johnston et al. 2005; 
Fellhauer et al. 2006; Law et al. 2009; Law & Majew- 
ski 2010). However it turned out that the conclusions 
of these studies depended on the specific set of obser- 
vational constraints considered, and it is now generally 
accepted the idea that none of the static-potential ax- 
isymmetric halo models considered is able to reproduce 
simultaneously all the available positional and kinematic 
data (see Yanny et al. 2009a; Law et al. 2009, for refer- 
ences and discussion) . In a recent contribution Law et al. 

(2009) anticipated that the adoption of triaxial halo mod- 
els can help to solve this problem: in Sect. 6.3 we briefly 
consider the N-body model they produced as a follow-up 
of that analysis (Law & Majewski 2010). In any case, it 
is quite clear that currently available models are far from 
perfect, and more detailed simulations are needed to ex- 
tract all the possible information on the Galactic DM 
halo from the Sgr Stream, as more (and more accurate) 
observational constraints become available. For exam- 
ple, Fellhauer et al. (2006) interpreted the bifurcation of 
the trailing arm giving rise to the A and B branches con- 
sidered here as produced by the precession between two 
subsequent orbits. As the implied amount of precession is 
relatively small, this, in turn, requires that the potential 
felt by Sgr should be nearly spherical. However the sim- 
ilarity between the two branches (in terms of distance, 
kinematics and stellar content) led Yanny et al. (2009a) 
to suggest that in fact the two branches are composed 
by stars lost at the same epoch, i.e. they are in the same 
orbital phase. In this case the separation between the 
two branches would not be related to orbital precession 
and would have nothing to say about the shape of the po- 
tential. In their recent analysis, Niederste-Ostholt et al. 

(2010) adopt the same view as Yanny et al. (2009a)^^. 
In the present contribution we provide very accurate 

distance estimates along the northern branches of the Sgr 
Stream as powerful constraints for future generations of 
Sgr disruption models that will include effects like halo 
triaxiality, dynamical friction, time-evolving Galactic po- 
tential, etc. In this section we discuss our flndings in 
comparison with the predictions of the three models by 
Law et al. (2005, L05 hereafter), just to show how power- 
ful accurate distance constraints can be in distinguishing 
between different models (some examples of such com- 
parisons have already been presented in Sect. 3.2 and 
Fig. 13). One of the main aims of the studies by L05 
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Fig. 20. — Comparison between the positions of the observed 
RC peaks and the predictions of Law et al.'s N-body model of 
the disruption of Sgr within an oblate DM halo in the X^^sgr vs. 
YQSgr plane (i.e. the approximate plane of Sgr orbit seen face- 
on). Heavier dots indicate model particles enclosed in our FoVs 
cones. To be consistent with Law et al. (2005) we converted our 
differential distances into absolute ones by assuming a true distance 
modulus of 16.90 for the main body of Sgr. 
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Fig. 21. — The same as Fig. 20, but for a N-body model within 
a spherical DM halo. 

and Johnston et al. (2005), was to use the existing ob- 
servations on the Sgr Stream to constrain the shape of 
the Galactic halo. For this reason the three models they 
provide^^ describe the final state of the evolution of a 
realistic progenitor of Sgr after a few orbits within a 
Galactic potential having a prolate, spherical, or oblate 
DM halo. For sake of simplicity, in the following we will 
refer to these models as to the Oblate (O), Spherical (S) 
and Prolate (P) models, respectively. In all the mod- 
els, each particle is fiagged according to the peri-galactic 
passage in which it become unbound from the main body 
of the galaxy. Here we refer to stars still bound or lost 



See also the discussion in Law & Majewski (2010). 
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Fig. 22. — The same as Fig. 20 and 21, but for a N-body model within a prolate DM halo. 



during the current peri-galactic passage as having p = 0; 
p = —1, —2, —3, —4 refers to particles lost one, two, three 
and four peri-galactic passages ago, respectively, p = Q 
stars are out of the range accessible with the fields con- 
sidered here, according to the LOS models. When we 
speak of "young" and "old" wraps of the Stream we re- 
fer to portions of the Stream whose population is domi- 
nated by particles lost in the most recent or less recent 
peri-galactic passages, respectively, on an age-scale en- 
compassing the last ~ 5 orbits, i.e. ~ 3 — 4 Gyr. 

In Figg. 20, 21 and 22 the three models are com- 
pared with the positions of the observed RC peaks in 
the Xq^sqt vs. 5^0, Sgr plane^^, as in Fig. 11, above. It is 
immediately apparent from Fig. 20 that the trend traced 
by our primary peaks rules out the O model, that fails 
to reproduce the most prominent branch of the Stream 

^^ For brevity, in the following we will drop the 0, Sgr index 
any time we found this convenient; XQ^Sgr, Y@,Sgr and X, Y are 
interchangeable. For the same reason, the values of X, Y must be 
always intended as expressed in kpc, even if not explicitly stated. 



seen in SDSS data, i.e. the portion of the leading arm 
descending from the North Galactic Pole (while there is 
some agreement for the - putative - nearby portion of the 
trailing arm). The case of the spherical model is similar, 
even if the disagreement between observations and model 
prediction is less severe (Fig. 21). 

The comparison with the prolate model is the most 
interesting and we take it also as the occasion to de- 
scribe the trends found in our data in a deeper detail. 
For a more fruitful discussion we provide Fig. 22 in a 
larger format with respect to its analogs for the oblate 
and spherical models (see also Fig. 13). It should be 
stressed that, in the following, we interpret the coherent 
structures we have detected using this specific model as 
a guideline. For an example of a different interpretation 
see Sect. 6.3, below. There are several features worth 
noticing in Fig. 22: 



1. For X, 



G,Sgr 



> 



-7 kpc the agreement between the 
positions of our primary branch-A peaks and the 
portion of the leading arm going from {X, Y) ~ 
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Fig. 23. — Comparison of observed peaks positions with Law et 
al.'s models in the XQ^Sgr vs. ZQ^Sgr plane (i.e. perpendicular 
to the orbital plane of Sgr, that is seen edge-on). The symbols are 
the same as in Fig. 20. 

(16,-40)19 to {X,Y) ~ (-7,-28) is excellent. 
Also primary branch-B peaks follow the same trend 
thus confirming that the two structures lie at the 
same distance (see also Fig. 13, and Belokurov et al. 
2006; Fellhauer et al. 2006; Yanny et al. 2009a). Ac- 
cording to the considered model, this part of the 
leading arm is dominated hy p = —1 particles up 
to XQ^sgr — 0? and by a mix of p = —2 and p = —3 
particles for Xq^sqt ^ 0. 

2. Several detections seem to extend the path of the 
arm down to (X,Y) = (-15,-15), possibly suggesting 
a slightly less elongated shape of the arm with re- 
spect to the model predictions. The coherence of 
the structure is less clear in this region: the models 
predict that various wraps cross here and this may 
be source of some confusion. 

3. The model predicts the presence of a more ancient 
(mostly populated hy p — —2 and p — —i particles) 
and wider wrap running nearly parallel to the por- 
tion of the leading arm described above, but behind 
it. This structure has been detected in branch-B, 
where one flag=2 point, at {X,Y) ~ (-5,-44), 
in coincidence with a branch-A detection, and a 
flag=3 detection at {X,Y) ~ (-14,-35)]. These 
points appear to trace the outer edge of this wrap 
, as depicted by the considered model. On the 
other hand, there is no detection (in any branch) 
for X > 0, i.e. where the detection of the second 
wrap should be easier, according to the model, as 
the separation from the inner wrap increases with 
X and the feature is denser and narrower in that re- 
gion. This lack of detection seems confirmed by the 
independent results of Belokurov et al. (2006), that, 
however, detect the most distant wrap at X < 
only in the direction of Branch A. To have a deeper 

^'^ In the following we drop the 0, Sgr indices, for brevity. The 
unity is always kpc. 



insight into this problem in Fig. 24 we provide a 
direct comparison between observations and model 
at the SCP level, as done in Fig. 12. Here we com- 
pare the observed SCP of the FlA, F2A and F3A 
fields with the SCPs obtained from the model in the 
considered los for particles lost one, two, three and 
four peri-galactics ago. From the upper-right panel 
it is clear that the dense X > part of the outer 
wrap, produced hy p = —2 particles in the model, 
has no counterpart in the observed SCPs and would 
be easily detected if actually there. On the other 
hand, the sum of the relics having p = -1,-3, 
and p — —4 provide a satisfactory match to all 
the observed peaks. This suggests that there is a 
real mismatch between the L05 P model predictions 
and our observations in this part of the halo. We 
note that the spherical model suffers from the same 
problem, while the oblate model does not predict a 
strong signal at that position, but it fails to match 
all the observations at X < for this wrap. 

A coherent series of detections lying at nearly con- 
stant y ~ 25 kpc, traced from X ~ 8 to X ~ —14 
in both branches, traces a filamentary structure 
that is identified here for the first time. Isolated 
detections with M giants and RR Lyrae were pre- 
viously reported at A ~ 295° (Majewski et al. 2003; 
Vivas & Zinn 2006; Ivezic et al. 2000a). This fea- 
ture matches quite well a wrap of the trailing arm 
that is present in all the L05 models; it can be ap- 
preciated from Fig. 13 and Fig. 22 that the agree- 
ment with the P model is very good. For X < — 14 
however the positions of the peaks do not trace the 
model prediction anymore. This apparent discon- 
tinuity along this branch cannot be (only) due to 
the distance effects discussed in Sect. 2.3 as the dis- 
tance is expected to increase and the sensitivity of 
the method should increase accordingly. Moreover 
we are able to detect peaks both more and less dis- 
tant than the position predicted by the model along 
these los. This feature has no counterpart in the 
triaxial halo model discussed in Sect. 6.3, below. It 
is clear that additional information is needed to un- 
derstand better the nature of this structure, from, 
for example, the kinematics of member stars. 

There are primary and tertiary branch-B detec- 
tions, plus one tertiary branch-A detections, trac- 
ing a feeble (but coherent) spur from an ancient 
(p = —3, —4) wrap, predicted by the model to 
arch between (X,F) ~ (-10,-15) and (X, F) ~ 
(—3, —22). As far as we know this is the first de- 
tection of this nearby portion of the Stream. A 
couple of primary branch-B detections (and a ter- 
tiary branch-A detection) may trace similar sub- 
structures on the near side of the constant-distance 
portion of the trailing arm (see Sect. 6.3 for an al- 
ternative interpretation) . 

There are a couple of other cases of slight distance 
mismatches between branch-A and branch-B de- 
tections, occurring, however in the region around 
(X, y) ^ (-10,-25) where different wraps of the 
Stream cross each other. It may be challeng- 
ing to disentangle the various contributions based 



24 



Correnti et al. 



on distances alone. A more interesting case is 
provided by the two pairs of detections around 
{X,Y) ^ (—23,-30), a region where the model 
predicts only feeble structures and branch-B detec- 
tions are clearly more nearby than branch-A ones. 
It is intriguing to note that the few particles of 
the model lying in this region are not uniformly 
distributed but appear to form two approximately 
parallel tiny bridges that reasonably reproduce the 
observed pattern. Also in this case this is the first 
detection of such structures. 

7. Both the P and S models by LOS predict the pres- 
ence of a fairly dense and narrow wrap composed by 
p = —3 and p — —4 particles crossing the accessi- 
ble range of the X,Y plane from (AT, Y) ~ (-7, -7) 
to {X,Y) ^ (5,-12), where it emerges from the 
d < 12 kpc zone of insensitivity of our method 
(the triaxial model briefly discussed in Sect. 6.3 
also displays a similar feature). Here we have a 
primary detection from the SCP of FIA, that cur- 
rently is the first detection of this nearby wrap of 
the leading arm^'^. A detailed exploration of the 
d < 12 kpc zone would require a different kind 
of analysis, hence it is postponed to a future con- 
tribution. We note, however, that Monaco et al. 
(2007) studied the chemical composition and the 
kinematics of a small sample of M giants that can 
be attributed to this nearby wrap. 

8. Two Branch-B detections, located at (A, Y) ^ 
(-3,-16) and {X,Y) - (+3,-18), F5B and F2B, 
respectively, do not seem to match any signifi- 
cant structure of the spherical and prolate LOS 
models; the primary one (that with positive X) 
is marginally consistent with the part of the lead- 
ing arm plunging toward the Sun of the oblate 
model (but see Sect. 6.3, below). As anticipated 
in Sect. 4.1, their position [(RA,Dec) = (195°, +16°) 
d - 18 kpc, and (RA,Dec)=(190°, +18.S°) d - 19. S 
kpc, respectively] is fully compatible with the outer 
fringes (i.e., the high Galactic latitude edge) of 
the nearby overdensity S297+63-20.S, discovered 
by Newberg et al. (2002) and discussed in detail 
in Newberg et al. (2007). It is unclear why we do 
not detect the structure in other adjacent fields, or 
in the corresponding branch A fields. This may 
be due to the intrinsic weakness of the RC sig- 
nal from these nearby features, or it may reflect 
a high degree of complexity of the sub-structures, 
as suggested in the analyses by Keller et al. (2009) 
and Vivas et al. (2008). Newberg et al. (2007) 
provided positional and kinematic evidence argu- 
ing against the association of S297+63-20.5 with 
the Sgr Stream, that was originally proposed by 
Martmez-Delgado et al. (2007) and cannot be com- 
pletely ruled out at the present stage (see also the 
discussion in Law & Majewski 2010). However 

^^ While the distinction between leading and trailing amis is easy 
and sensible for particles lost in the latest two peri-galactic pas- 
sages, it becomes increasingly blurred for Stream wraps dominated 
by more ancient relics, as a particle can reach the same position 
in these parts of the Stream both from the leading and from the 
trailing sides of the tidal tails. 



Fig. 22 provides further support for the conclu- 
sions by Newberg et al. (2007): the peaks detected 
here do not present any continuity with the main 
branches of the leading and trailing arms of the 
Stream as traced in the present analysis (but see 
also Sect. 6.3, below). Our data suggest that the 
leading arm crosses the Galactic plane at ~ 10 kpc 
from the Sun, toward the Anticenter, in agree- 
ment with Newberg et al. (2007) and Seabroke 
et al. (2008). On the other hand, the identifica- 
tion of S297+63-20.5 with the Virgo Stellar Stream 
(VSS, Duffau et al. 2006; Vivas et al. 2008) seems 
likely, while the relationship between VSS and the 
Virgo Over Density (VOD, Juric et al. 2008; New- 
berg et al. 2007; Keller et al. 2009) is less certain 
(see Newberg et al. 2007; Keller 2009b). We are 
currently following up these possible detections of 
S297+63-20.5/VSS in FSB and F6B (also looking 
for the structure at lower latitudes). If confirmed, 
they would provide the first detection of RC stars 
in these structures, in analogy with the cases of 
Boo HI discussed in Correnti et al. (2009). RC 
stars may provide new insights on the nature of 
complex series of structures recently identified in 
the direction of Virgo (Keller et al. 2009; Keller 
2009b). 

All the features and correlations with the P model de- 
scribed above can be seen even more clearly and directly 
in Fig. 13, that provides the most natural way to com- 
pare our measures with models. For example, the match 
between two weak model structures described at point 
6, above, and our detections can be very clearly appre- 
ciated in that plot, at 230° < A < 24S° and d ~ 37 kpc. 
The linear trend of increasing distance with decreasing A 
of the two parallel sets of observed points is very nicely 
matched by corresponding filaments of particles in the 
model. 

Fig. 23 shows that the overall morphology of the three 
models is remarkably similar in the Xq^sqt vs Zq^sqt 
plane and reproduces the general trends of the data (ex- 
cept for the Oblate model, that predicts a total lack of 
particles for branch-B detections at Xq^sqt > 0, at odds 
with observations). A more detailed analysis is beyond 
the scope of the present paper. On the other hand, we 
must conclude, from the results summarized above, that 
the prolate model by Law et al. (200S) is the one (among 
those considered here) providing the best match to the 
positional data considered here. It should be stressed 
that with this we do not intend to say that a prolate halo 
model is favored by our data, as the comparison was lim- 
ited to just three very specific models that are already 
known not to be able to fit all the positional and kine- 
matical observational constraints available (Law et al. 
2009). In particular it should be recalled that the avail- 
able radial velocities of Stream stars seems to favor pro- 
late models (Helmi 2004), while the angular precession 
of the leading arm with respect to the trailing arm favor 
spherical or slightly oblate models (LOS; Johnston et al. 
200S; Newberg et al. 2007; Prior et al. 2009a, and ref- 
erences therein). We simply note that any future model 
intended to fit all the observed characteristics of the Sgr 
Stream must have a spatial structure very similar to that 
of the prolate model by Law et al. (200S), at least in the 
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Fig. 24. — Comparison of the observed SCPs in the field FIA, F2A, F3A with the synthetic SCPs obtained by adding the (arbitrarily 
normalized) los histogram from the prolate N-body model (as in Fig. 12 above) to the best-fit background model of the considered field. 
Each of the four (triple) panels report the synthetic SCPs including only N-Body debris stripped from the main body one, twro, three or 
four orbits ago, going from the upper left, to the upper right, to the lower left and to the lower right panels, respectively. 

actual width. This expected effect is clearly confirmed 
in Fig. 25, where observed FWHM are always equal or 
larger than their model counterparts. The disentangle- 
ment of overlapping structure may be also problematic, 
as it is unavoidably performed in different ways in the 
observed SCPs and in the N-body models. 

To minimize the possible ambiguities associated with 
the collapse of complex structures along the los into a 
single FWHM measure (see Sect. 3), especially in regions 
where different wraps cross one another, in Fig. 25 we 
limit our comparison to X > —10 kpc peaks tracing the 
two main wraps (leading and trailing arms) that are 30 
kpc apart at Asgr — 290° and cross each other at Asgr — 
265°, and we consider only primary and secondary peaks. 

The most interesting and sensible comparison is be- 
tween the trends of the FWHM as a function of orbital 



portion of space sampled by our study, unless an alter- 
native origin is assumed for the d ~ 25 kpc structure we 
tentatively interpreted as the trailing arm. 

6.1. Trends of depth as a function of orbital azimuth 

In line with the above discussion, in Fig. 25 we com- 
pare the FWHM along the los described in Sect. 3.5 with 
those measured from the distribution of particles of the 
Prolate model of L05, along the same los . The following 
discussion is mainly intended to illustrate the possible 
use of the derived FWHM. It should be considered that 
there are additional sources of uncertainty affecting this 
comparison, associated with the the measure of FWHM 
in models. For example, the measured width depend on 
the actual number of particles of the model, the limited 
number of particle may lead to underestimates of the 
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Fig. 25. — Comparison between the trends with Agg,, of the 
observed FWHM along the los (the same symbols as Figg. 20, 21 
and 22, above, and the predictions of the Prolate model by LOS 
(open triangles). The comparison is limited to the los from Fl to 
F7 (i.e., those providing the cleanest tracing of the leading and 
trailing arms). The upper row of panels refers to detections in the 
leading arm, the lower row to detections in the trailing arm. The 
left and right columns refer to branch A and branch B, respectively. 
The field numbers increase (from Fl to F7) from right to left, as 
in Fig. 13 and Figg. 20, 21, 22. 

azimuth. The upper panels of Fig. 25 show that the 
observed and predicted trends for the leading arm are 
indeed similar, both in direction and in amplitude, for 
both branches. The agreement of the absolute values of 
the FWHM is also satisfying (within a factor of ^ 2), 
with four Branch A and one Branch B detections closely 
matching the model predictions. The FWHM of the con- 
sidered detections from FlA and F3A give some reason 
of concern, as they break the continuity of the observed 
trend: this may suggest that there may be some unre- 
solved structure in these peaks. Alternatively we have 
to accept variations of a factor of ~ 2 as due to the un- 
certainty inherent to the adopted method of estimating 
FWHM. The overall agreement is reasonable also for the 
putative trailing arm. 

It is interesting to note that the different trends ob- 
served in the two branches of the leading arm are repro- 
duced by the P model that do not present any bifurcation 
(see also Fig. 26, for a similar behavior in the A vs. den- 
sity trend in the leading arm). 

6.2. Trends of density as a function of orbital azimuth 

In strict analogy with the analysis described in the 
previous subsection, in Fig. 26 we present the comparison 
of the observed and predicted trends of the stellar density 
(see Sect. 3.4) as a function of A. The measured density 
is compared with the density of particles in the same 
wrap of the P model. The density scale of the model has 
been multiplied by the arbitrary factor 2.5, to achieve a 
reasonable normalization with the observed values. As 
in Sect. 6.1 the comparison presented is just intended 
as illustrative of the possible use of these numbers, and 
it is limited to the cleanest portions of the leading and 
trailing arms, at AT > —10 kpc. 

For the leading arm the match between the overall ob- 
servations and the models is acceptable, in particular for 
branch A. The highly discrepant point at A ~ 263° is 
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Fig. 26. — The same as Fig. 25 for the density of RC stars as a 
function of A. An arbitrary normalization factor of 2 is applied to 
all the density values from the theoretical model. 

associated with an especially complex SCP model, with 
three overlapping peaks (F5B): for this reason we are in- 
clined to ascribe the discrepancy to an erroneous density 
estimate. 

On the other hand, while the model predicts low or 
even negative gradients of density with increasing A, the 
observations show a very strong positive gradient, sim- 
ilar in both branches^^. This is an obvious example of 
the kind of constraints that can be achieved with these 
data: in principle, any fully successful model of the dis- 
ruption of Sgr must also reproduce a density gradient 
similar to the observed one. However it has to be taken 
into account that the available models are intended to 
describe the dark matter halo in which the baryonic part 
of the galaxy is embedded. While, for example, stars and 
DM particles in the Stream should not greatly differ in 
their kinematical and positional properties, their density 
would follow the same trends only if mass strictly follows 
light also in tidal tails, which is very unlikely to be the 
case (see Peharrubia et al. 2008, and references therein). 

6.3. The triaxial model by Law & Majewski 2010 

When the present manuscript was ready for submis- 
sion, a preprint was posted (Law & Majewski 2010), fol- 
lowing up the preliminary analysis by Law et al. (2009). 
In that study a new N-body model of the disruption of 
Sgr within a triaxial Galactic potential, is shown to pro- 
vide a reasonable match to most of the existing observa- 
tional constraints. In particular, the new model repro- 
duces the distance vs. A trend reported by Bel06 for the 
main wrap of the leading arm, the precession between 
the leading and trailing arms, and it matches the exist- 
ing sets of kinematic measures. 

As we have stressed before, it is far beyond the scope 
of the present analysis to find out which is the best avail- 
able model. However it is worth showing the comparison 
between our results and this new model, for complete- 

■^^ It has to be recalled that the bifurcation in two branches is 
an observed property of the main wrap of the leading arm. There 
is no reason to discuss other wraps as divided into two branches: 
here this is merely incidental, due to adopted distribution of the 
observed fields that were chosen to trace the bifurcation. 
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Fig. 27. — Comparison with the new "triaxial halo" (T) model 
by Law & Majewski (2010): distance vs. orbital azimuth. Symbols 
are the same as in Figg. 20, 21 and 22. Our measures have been 
rescaled to match the distance scale adopted by Law & Majewski 
(2010), that implies a distance to the main body of Sgr d = 28 kpc 
(i.e. (m — M)o = 17.23, instead of (m — M)o = 16.90, as adopted 
by L05). 



Triaxial Halo Model 
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Fig. 28. — Comparison with the new "triaxial halo" (T) model 
by Law &; Majewski (2010): the Xq^Sqt vs. YQ^Sgr plane, with 
the same assumptions and symbols as Fig. 6.3. 

ness and (above all) as a very instructive example of how 
the interpretation of observed features may depend on 
the considered model (see Sect. 7 for discussion). From 
the inspection of Fig. 27 and Fig. 28 it can be concluded 
that the triaxial halo model provides a reasonable match 
to distance gradient of the main wrap of the leading arm, 
over the whole range of A covered by our data. This is not 
unexpected as the model is found to fit the observations 
by Bel06 in this regard, and we are in good agreement 
with these authors. The same is true for the sparse de- 
tections behind the main wrap, that were also found by 
Bel06. It is interesting to note that, in a similar fashion 
to the P and S models, the T model predicts a remark- 



able increase of the density of this wrap for A > 275° 
{X > 0) that is not observed, neither in the present 
work or by Bel06. Moreover, the T model does not seem 
to display the narrow and dense structure of the main 
wrap of the leading arm that in the P model appears to 
match so well our coherent set of primary detections in 
that region. We postpone a detailed comparison between 
the observed structure along the los and the predictions 
of the P and T model to a future contribution: here we 
limit the discussion to the main features of the models 
(i.e. trends of distance with orbital azimuth). 

The new model makes predictions very similar to those 
of the P and O models also regarding the nearest wrap of 
the Stream, running across the whole range of A sampled 
by our data. However it should be noted that it predicts 
this wrap to lie below our sensitivity limit at any A, in the 
observed range, thus it is unable to match the observed 
points at A ~ 287° and d ~ 13 kpc, at odds with the 
P model. The T model presents a very coherent narrow 
wrap of the trailing arm running at nearly constant d ~ 
20 kpc from A ~ 290° to A ~ 250°, then it begins to 
bend gently toward d ~ 25 kpc from A ~ 235° where 
it crosses with the leading arm. This feature matches 
very nicely the nearby (d < 25 kpc) detections that we 
tentatively attributed to S297-H63-20.5/VSS and to an 
ancient spur of the leading arm, in the comparison with 
the P model described in detail in Sect. 6, above. 

On the other hand, the coherent structure we detect 
from A ~ 260° to A ~ 250°, that we interpreted as 
a wrap of the trailing arm, is not present in the T 
model. The same is true for the d > 40 kpc structures at 
A < 245°. The orbital path of the simulated Sgr galaxy 
matches also these structures, so it is not excluded that 
they may correspond to very ancient wraps. However, it 
has to be noted that the T model is the remnant of the 
evolution of a Sgr progenitor for ~ 8 orbits (not just ~ 4 
as for S,0, and P) models, thus it should include wraps 
populated from more ancient stripping events than the 
S, O, and P models. 

In conclusion, while the P model still appears to pro- 
vide a more thorough match of the observed structures, 
the T model provides a promising alternative that de- 
serves to be investigated in further detail. Not surpris- 
ingly, the mere comparison with our own (limited) data- 
set shows that both models need to be refined. 

7. SUMMARY AND CONCLUSIONS 

We have used RC stars to trace the long tidal tail of 
the Sgr dSph galaxy in the portion of the Northern sky 
sampled by the SDSS-DR6. Structures along the line-of- 
sight are identified as peaks in the (otherwise smoothly 
increasing) /- and y-band SCPs of color-selected samples 
of candidate RC stars, from ~ 5° x 5° fields covering the 
whole extension of the two main branches (A and B) of 
the Sgr Stream identified by Belokurov et al. (2006) in 
the same dataset. Any other part of the Stream in addi- 
tion to these branches is expected to lie (approximately) 
in the same plane, i.e. it should be visible in the consid- 
ered fields. The analysis was focused on obtaining the 
most accurate and reliable distances to all the wraps of 
the Stream that we were able to detect. 

Many significant peaks were consistently found in both 
the SCPs of several fields. The observed SCPs were 
modeled as a series of Gauss curves (one for each peak) 



28 



Correnti et al. 



superposed to a polynomial accounting for the smooth 
fore/background population. For each significant peak 
we derived a purely differential estimate of its distance 
(with uncertainties < 10%), an estimate of the FWHM 
along the los , and an estimate of the associated density 
of RC stars attributable to the considered structure. AU 
the derived quantities are provided in Tab. 2 as powerful 
constraints for the new generations of models of the dis- 
ruption of the progenitor of Sgr dSph within the Milky 
Way halo. 

To illustrate the potential of our measures in that con- 
text we compared them with the three models made pub- 
licly available by LOS. These provide a realistic realiza- 
tion of the present epoch configuration of particles that 
were originally bound to a progenitor similar to Sgr that 
was evolved for ~ 4 orbital periods within a static Galac- 
tic potential with different degrees of flattening (a spher- 
ical, oblate and prolate halo, respectively). The models 
(and in particular the Prolate halo one, that matches 
well most of our observations) are also used as guidelines 
for the interpretation of our results. The great complex- 
ity of a structure like the Sgr Stream, multiply-wrapped 
around the Galaxy, requires a process of convergence be- 
tween models and observations: the latter must constrain 
models but the former are indispensable to re-conduce 
such a complexity to a single structure (see Sect. 6.3). 

Our technique resulted in higher- accuracy distance es- 
timates with respect to previous studies, and demon- 
strated high sensitivity to feeble structures. However, 
the sensitivity is easily destroyed by contamination from 
Galactic sources: for these reasons we had to limit our 
survey to 6 > 50°, while other (more abundant) tracers 
are able to follow the Stream down to 6 > 30° . The over- 
all agreement with previous analyses is good (see Sect. 5). 
Finally, and most importantly, our method proved es- 
pecially efficient in the detection of (relatively) nearby 
structures. In the following we summarize and briefiy 
discuss the main conclusions of the present study, taking 
Fig. 13 and Fig. 22, as references. 

• For A > 255° (X > -10 kpc) the leading arm of the 
Stream is cleanly and coherently detected in both 
branches, going from d = 43 kpc at A ~ 290° to d = 
30 kpc at A ~ 255°. This is in full agreement with 
the results obtained with other tracers (Newberg 
et al. 2007; Niedcrste-Osthoh et al. 2010). This 
portion of the leading arm is the most unambiguous 
and robustly constrained. 

• In the same range of A (and X) a remarkably coher- 
ent structure is also very clearly detected at nearly 
constant distance from us, d ~ 25 kpc. According 
to the S and P models by L05 this can be inter- 
preted as a wrap of the trailing arm, while it has 
no obvious counterpart in the recently presented 
T model (Law & Majcwski 2010). The P model 
matches the observed structure very well. Previous 
detections of this wrap were reported only around 
A — 295° (see Majewski et al. 2003, for discussion 
and references). 

• The comparison with the L05 models strongly sug- 
gest that the run of the relative distance as a func- 
tion of A of the two wraps described above has 



a strong power in discriminating between differ- 
ent models of the Stream. In particular the S and 
O models by L05 clearly fail to reproduce the ob- 
served trends. On the other hand the P model 
reproduces the trend nearly perfectly. 

• Weak detections of a further, more distant wrap 
(running parallel to the leading arm, in the same 
range of A as above) were also obtained. These 
support similar results by Belokurov et al. (2006). 
An enhancement of the density of this wrap at A > 
275°, predicted by the S, P and T models, seems to 
be excluded by the present analysis (in agreement 
with Bel06). 

• Turning to the A < 255° (X < -10 kpc) portion of 
the survey, this is characterized by a very complex 
structure, partly due to the crossing of multiple 
wraps predicted to occur in this region by all the 
models. Hence the interpretation of these struc- 
tures is less straightforward, and must be consid- 
ered as tentative. However, the P model appears 
to provide a reasonable match to all the detections 
in this region: for these reasons we adopt it as a 
guideline for our best-effort interpretation of the 
data (see Sect. 6.3 for an alternative view). 

• The leading arm seems to be traced beyond A ~ 
255°, continuing its trend of linear decrease of its 
distance down to d ~ 20 kpc at A ~ 220°. Extrap- 
olating from the observed trend one would expect 
the arm to cross the Galactic disk at ~ 10 kpc 
from the Sun, in agreement with the conclusions 
by Newberg et al. (2007); Seabroke et al. (2008); 
Law & Majewski (2010). The degree of coherence 
of the detections in this portion of the leading arm 
is lower, suggesting the possible presence of further 
(unresolved) substructure or due to higher uncer- 
tainties associated to weaker and overlapped struc- 
tures. 

• In the same region, the continuation of the trailing 
arm is coherently traced where predicted by the 
P model up to A > 240° {X > -15 kpc). For 
A < 235°, in particular, we lack any detection cor- 
responding to the well defined structure predicted 
by the model (the same is true for the T model). 
On the other hand, coherent detections are ob- 
tained behind the main wrap of the trailing arm 
as predicted by the P model for A < 240°. These 
detections may indicate a different shape for that 
portion of the trailing arm. However, as discussed 
above, they match two more feeble structures run- 
ning parallel to the main arm. It is obvious that the 
P model is not adequate to fit all our observations, 
in spite of the good overall match. 

• The most nearby detections are the more difficult 
to interpret robustly. However the single primary 
detection at d ~ 13 kpc and A ~ 287° (just beyond 
the d < 12 kpc "zone of avoidance" of our tech- 
nique) matches the prediction of all the three L05 
models, as well as for the model by Law & Majew- 
ski (2010). For this reason we are quite confident 
to have detected for the first time the nearest wrap 
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of the leading arm. We are currently following up 
this finding, to check if the predicted d ^ 10 kpc 
wrap can be detected also in other los . 

• The three detections at d ^ 20 kpc and A ~ 245° 
are matched by a spur of the P model. The two de- 
tections at d ~ 18 kpc and A > 260° have been ten- 
tatively ascribed to the S297+63-20.5/VSS over- 
density. The T model matches very well all of these 
detections with a single narrow wrap of the Sgr 
trailing arm. However Law & Majewski (2010) con- 
firm that the kinematics predicted by their model 
toward VSS is markedly different from what ob- 
served by Duffau et al. (2006) and Newberg et al. 
(2007). 

• The overall trends of FWHM along the los as a 
function of A of the P model provide a reason- 
able match to our primary detections of the leading 
arm. It is especially interesting to note that the 
model reproduces the different trends encountered 
in the two branches, even if it does not produce the 
observed bifurcation. This seems to provide fur- 
ther support to the view (adopted by Yanny et al. 
2009b; Niederste-Osthoh et al. 2010) that branch A 
and branch B are substructures within the same 
wrap of the Stream, and not different wraps as pro- 
posed by Fellhauer et al. (2006). However, a lim- 
ited set of tests performed on intra-branches fields 
suggests that the bifurcation in the Dec direction 
shown by Belokurov et al. (2006) is real. Proba- 
bly a deeper, thorough and independent analysis 
of the Dec structure of this wrap of the Stream is 
warranted (see Law & Majewski (2010) for possible 
alternative explanations). 

• The observed trends of density as a function of 
RA along the leading arm (branches A and B) are 
in fair agreement with those by Niederste-Ostholt 
et al. (2010). Our estimates of the total luminosity 



per kpc at any given RA are lower than theirs by 
a factor of ^ 4 — 5. 

• Kinematic follow up of the newly identified struc- 
tures is clearly urgent. Carrell & Wilhelm (2010) 
recently demonstrated that this can be carried on 
using exactly the same tracer stars, i.e. RC stars. 
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